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1.  SCOPE. 

This  Test  Operating  Procedure  (TOP)  presents  considerations  and  techniques  for  developing 
Laboratory  Vibration  Test  Schedules  (LVTS)  that  can  be  utilized  to  simulate  field  vibration 
environments  on  a  vibration  table.  Laboratory  vibration  tests  are  used  extensively  in  lieu  of 
more  time-consuming  and  less  cost  effective  field  exposure  tests.  This  TOP  specifically 
addresses  random  vibration  testing  controlled  to  frequency-domain  vibration  spectra  and  is 
intended  to  address  multiple  “exciter”  (also  referred  to  as  “shaker”  or  “actuator”)  scenarios  with 
the  emphasis  on  mechanical  multiple  degree-of-freedom  (MDOF)  scenarios.  There  is  a 
significant  increase  in  complexity  between  single-exciter/single-axis  (SESA)  and  multiple- 
exciter/multiple-axis  (MEMA)  testing  in  terms  of  both  mechanics  and  control.  MEMA  specific 
issues  ranging  from  definitions  and  nomenclature  consistency,  to  data  analysis  techniques,  will 
be  addressed. 

2.  FACILITIES  AND  INSTRUMENTATION. 

2.1  Facilities. 

The  development  of  a  LVTS  will  require  access  to  the  test  item  of  interest  (or  a  dynamically 
equivalent  surrogate),  access  to  the  carrier  vehicle,  appropriately  placed  transducers,  signal 
conditioning  and  data  acquisition  hardware,  and  a  controlled  environment  for  collecting  input 
data  (e.g.,  a  road  course  for  wheeled  and/or  tracked  vehicles,  waterway  for  watercraft,  airspace 
for  aircraft,  rotorcraft,  and/or  spacecraft). 

2.2  Instrumentation. 

a.  LVTSs  are  generally  defined  in  terms  of  acceleration  units.  The  transducer  of  choice 
for  making  acceleration  measurements  is  an  accelerometer.  This  TOP  will  address  LVTS 
development  in  terms  of  acceleration. 

b.  It  is  strongly  recommended  that  the  same  model  of  accelerometer  and  signal 
conditioning  is  employed  at  all  instrumented  locations  to  preserve  phase  characteristics  during 
both  the  field  acquisition  and  laboratory  test  phase  of  any  MDOF  test.  Refer  to  the  guidelines  in 
Military  Standard  (MIL-STD)-810G'  and  Institute  of  Environmental  Sciences  and  Technology 
(IEST)  Recommended  Practice  IEST-RP-DTE012.22  for  recommended  accuracy  of  the 
transducers  and  associated  signal  conditioning. 

3.  REQUIRED  TEST  CONDITIONS. 

The  primary  function  of  Vibration  Schedule  Development  (VSD)  is  to  combine  vibration 
measurements  of  numerous  events  that  collectively  represent  an  item’s  lifetime  vibration 
exposure  (or  some  predefined  subset  thereof)  into  a  manageable  set  of  LVTS  representing  the 
equivalent  exposure.  The  most  dynamically  accurate  method  to  reproduce  the  full  exposure 
would  be  to  sequentially  vibrate  the  system  to  all  the  individual,  uncompressed  events 

*Superscript  numbers  correspond  to  those  in  Appendix  E,  References. 
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representing  its  full  lifecycle.  However,  such  an  approach  is  generally  not  feasible  from  both 
schedule  and  economic  perspectives  and  some  compromises  must  be  made  to  realize  the  benefits 
of  testing  in  the  laboratory.  Time  compression  techniques  based  on  fatigue  equivalency  are 
typically  employed  such  that  vibration  testing  can  be  performed  in  a  timely  and  economic 
manner.  North  Atlantic  Treaty  Organization  (NATO)  Allied  Environmental  Conditions  Test 

•5 

Publication  (AECTP)  240,  Leaflet  2410  ,  provides  general  guidance  for  developing  accurate 
representations,  and  issues  that  should  be  considered  during  the  VSD  process  for  the  SESA 
scenario.  This  TOP  expands  upon  the  discussion  in  Leaflet  2410  to  address  the  general  multiple 
exciter  test  scenario.  Discussions  will  be  limited  to  random  LVTS  development.  At  the  time  of 
this  publication,  no  commercially  available  multiple-input  multiple-output  (MIMO)  solutions 
exist  for  swept  narrowband  random  on  random  (NBROR)  or  sine-on-random  (SOR).  Refer  to 
MIL-STD-810G,  Method  525,  Procedure  I,  for  a  Time  Waveform  Replication  (TWR)  based 
MDOF  solution  for  swept  NBROR  and  SOR  requirements. 

3.1  Test  Configurations. 

The  MIMO  random  vibration  test  problem  can  refer  to  several  configurations.  One  configuration 
is  multiple  exciters  driving  a  single  test  item  in  one  axis.  This  configuration  is  often  used  for 
large  test  items  too  large  for  a  single  exciter.  A  second  configuration  is  the  excitation  of  a  single 
test  item  with  multiple  exciters  in  more  than  one  axis.  Linear  displacements  along  defined 
directions  are  referred  to  as  translation  degree-of-freedom  (DOF)  and  angular  displacements 
along  those  same  directions  are  referred  to  as  rotation  DOFs.  Up  to  six  DOFs  exist  for  a  rigid 
body  (i.e.,  X-,  Y-,  Z-translations  and  roll,  pitch,  yaw  rotations).  In  some  cases,  additional  DOFs 
can  be  excited  due  to  elastic  deformations  of  the  test  article. 

3.1.1  Basic  Representation  of  a  MIMO  System. 

All  MIMO  test  systems  are  discussed  using  a  common  description  in  tenns  of  matrix 
equations2'4'5.  A  simplified  version  of  the  general  MIMO  random  vibration  test  problem  can  be 
generalized  in  Figure  1 .  The  complete  mechanical  system  is  characterized  by  the  power 
amplifiers  and  a  system  of  several  exciters,  on  which  is  mounted  a  single  test  article.  The 
response  of  the  test  article  is  monitored  by  a  vector  of  response  channels  (represented  as  {c}). 
Each  element  in  the  vector  is  typically  the  acceleration  time  history  from  a  single  accelerometer. 
In  theory,  other  types  of  sensors  could  be  used.  The  power  amplifiers  are  driven  by  a  vector  of 
electrical  drives  (represented  as  {d}),  generated  by  a  control  system.  Each  element  in  the  vector 
is  a  time  history  driving  a  single  shaker.  The  control  system  monitors  the  response  of  the  test 
item  {c},  and  attempts  to  produce  drive  signals  {d},  such  that  the  statistics  of  the  control  signals 
meet  some  criteria  as  specified  in  the  test  specifications. 
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Figure  1.  Basic  representation  of  a  MDOF  system. 

3.1.2  Generalized  Representation  of  a  MIMO  System. 

A  more  generalized  MIMO  system  is  shown  in  Figure  2.  A  system  under  test  is  driven  by  Ns 
shakers  resulting  in  the  response  of  Na  control  accelerometers.  The  accelerometer  data  are 
typically  structured  in  blocks.  Each  of  the  acceleration  records  will  then  be  a  vector  of  time 
samples.  Some  control  systems  then  provide  for  a  transfonnation  matrix,  Ta,  to  convert  the 
block  of  Na  accelerometer  time  histories  to  Nc  control  variables.  The  Spectral  Density  Matrix 
(SDM)  of  the  control  variables  is  then  estimated  from  the  current  block  of  data  and  previous 
data.  The  transformation  matrix,  Ta,  is  typically  a  constant  independent  of  frequency.  In  theory 
the  transfonnation  matrix  could  be  applied  before  or  after  the  estimation  of  the  control  SDM. 

The  estimated  control  SDM,  C  is  then  compared  with  the  reference  SDM,  R,  and  a  conection  is 
computed  for  the  drive  SDM,  D.  The  drive  time  histories  {d}  are  then  computed  from  the  drive 
SDM,  D,  using  time  domain  randomization.  A  second  transformation  matrix,  Ts,  is  employed  to 
transform  the  Nd  drive  variables  into  Ns  shaker  drive  signals.  In  theory,  Ts  could  be  implemented 
before  or  after  the  transfonnation  into  the  time  domain.  One  advantage  of  placing  the 
transformation  in  the  frequency  domain  section  of  the  control  algorithm  is  that  the  matrix  could 
then  be  made  a  function  of  frequency.  Having  the  transfonnation  matrix,  Ts,  a  constant  assumes 
the  shakers  are  matched  and  the  desired  transfonnation  can  be  deduced. 
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Figure  2.  Generalized  representation  of  a  MDOF  system. 

3.2  Generalized  MDOF  Vibration  Control  Discussion. 

a.  A  general  discussion  of  the  MDOF  control  process  is  provided  for  insight  as  to  how 
the  MDOF  LVTS  will  serve  as  the  reference  in  the  control  process.  The  correction  can  be 
computed  in  several  ways.  One  method  is  to  compute  the  drive  from: 

D  =  ZRZ' 

where  the  system  impedance  matrix,  Z,  is  updated  as  new  infonnation  is  gathered,  or  a  modified 
reference  spectrum,  R ,  is  computed  based  on  the  error  in  the  return  spectrum.  The  initial  drive 
vector  is  typically  computed  using  the  above  equation  and  the  reference  SDM.  A  drive  signal 
error  can  also  be  computed  from: 


Dc  =  Z(R-C)Z' 

Sometimes  an  adaptive  correction  is  used.  Sometimes  a  combination  of  all  methods  is  used. 
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b.  The  transformation  matrices  are  often  called  the  input  and  output  transformation 
matrices.  One  should  be  careful  with  this  nomenclature  because  of  the  confusion  between  input 
and  output.  The  input  to  the  system  under  test  (voltages  to  the  power  amplifiers)  is  the  output  of 
the  control  system.  The  output  of  the  system  under  test  (the  accelerometer  measurements)  is  the 
input  to  the  control  system.  Paragraphs  4.4. 1  and  4.4.2  provide  the  nomenclature  employed  for 
input  and  output  transformations,  as  they  are  applied  within  this  document. 

c.  Minor  errors  in  the  matching  of  shakers  can  be  corrected  by  the  control  algorithm,  but 
major  mismatches  could  be  problematic.  The  time  domain  drive  signals  (represented  by  {s}), 
are  sent  to  the  shakers  completing  the  control  loop. 

d.  If  Ta  is  not  available,  then  Na  =  Nc  and  {a}  =  {c} .  If  Ts  is  not  available,  then  Nd  =  Ns 

and  {d}  =  {s}.  If  Nd  =  Nc,  the  number  of  control  variables  and  the  number  of  drive  variables  are 
the  same.  This  is  referred  to  as  square  control.  Square  control  is  the  most  common  control 
method.  If  Ns  >  Na  the  system  is  over-actuated  and  least  squares  approach  using  a  pseudo 
inverse  (pinv)  is  typically  used  to  determine  the  drive  signals.  If  Ns  <  Na  the  system  is  under¬ 
actuated  and  exact  control  of  the  control  SDM  is  often  not  possible.  In  such  cases,  some  kind  of 
average  control  is  usually  implemented.  Often  when  Ns  ^  Na  some  combination  of  the 
transformation  matrices  are  often  used  to  force  square  control,  Nd  =  Nc. 

e.  The  entire  mechanical  system  can  be  characterized  by  a  matrix  of  frequency  response 
functions  [H].  For  the  typical  case,  these  frequency  response  functions  will  have  units  of  g/V 
(acceleration  in  gravitational  units/volts  of  drive).  For  the  typical  case,  the  control  signals  are 
characterized  by  a  SDM.  The  diagonal  elements  are  the  autospectral  density  (ASD  or  PSD)  of 
the  control  signals.  The  off  diagonal  elements  are  the  cross  spectral  densities  (CSD)  between 
pairs  of  control  signals.  The  input  to  the  system  is  characterized  by  the  SDM  of  the  voltage  drive 
signals.  The  fundamental  relationship  between  the  drives  and  the  control  signals  is  given  by: 

C  =  HDH ' 

f.  The  complex  conjugate  transpose  is  denoted  by  [  ]’.  All  of  the  matrices  in  the 
equation  are  complex  functions  of  frequency.  The  spectral  density  matrix  is  Hennitian6,  i.e. 

Dy  =  Djt  *  where  ZT  *  is  the  complex  conjugate  of  D ,,  ,  and  Djt  is  an  element  from  a  spectral 

density  matrix.  Note  that  this  requirement  demands  that  the  diagonal  elements  are  real.  Note 
that  C  and  D  are  square  matrices;  they  have  the  same  number  of  rows  and  columns.  C  and  D  are 
the  same  size  only  if  H  is  square,  i.e.  the  same  number  of  inputs  and  outputs.  To  be  physically 
realizable,  the  SDM  must  also  be  positive  semi-definite.  This  requirement  will  be  discussed  in 
paragraph  4.5.2. 

g.  The  drive  spectral  density  matrix  is  converted  into  the  drive  time  histories  using  the 
method  of  time  domain  randomization4.  The  spectral  density  matrix  is  typically  estimated  using 
Welch’s  method7. 
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4.  TEST  PROCEDURES. 

VSD  requires  a  thorough  knowledge  of  the  dynamic  environment  to  which  the  test  hardware  will 
be  exposed  when  fielded.  This  knowledge  must  include  characterization  of  the  exposure  levels 
and  durations  for  all  relevant  conditions. 

4.1  Development  of  Mission  or  Lifetime  Scenario. 

The  duration  of  the  vibration  environments  can  be  derived  from  the  item’s  Life  Cycle 
Environment  Profile  (LCEP).  The  life  cycle  will  include  many  different  types  of  induced 
mechanical  environments  which  may  occur  while  the  materiel  is  being  handled,  transported, 
deployed  and  operated.  Although  all  the  induced  mechanical  environments  are  not  critical  in 
terms  of  generating  potential  damaging  response  amplitudes,  they  contribute  in  varying  degrees 
to  the  materiel’s  fatigue  damage.  All  expected  exposure  conditions  should  be  tabulated,  along 
with  corresponding  durations,  to  fonn  the  items  lifetime  “scenario”.  The  scenario  is  a  key 
parameter  in  the  development  of  any  vibration  schedule. 

4.2  Limitations. 

The  mechanical  degrees  of  freedom  (DOFs)  for  which  a  VSD  effort  is  capable  of  addressing,  is  a 
function  of  the  number  and  placement  of  the  transducers  employed  in  the  field  data  acquisition 
phase.  Similarly,  the  maximum  number  of  mechanical  DOFs  possible  to  reproduce  in  the 
laboratory  environment  is  a  function  of  the  number  and  placement  of  actuators  and  coupling 
hardware.  This  TOP  will  consider  the  general  case  for  VSD  development  in  which  the  reference 
SDM  will  be  defined  in  terms  of  the  six  classical  (3-translational  and  3-rotational)  rigid  body 
mechanical  DOFs.  In  the  event  less  than  six  mechanical  DOFs  are  being  considered,  the 
generalized  theory  is  easily  configured  to  address  the  motion  of  interest. 

4.3  Field  Data  Acquisition. 

When  in-service  measurement  data  have  been  obtained,  it  is  assumed  that  the  data  is  processed  in 
accordance  with  good  data  analysis  procedures,  as  in  Multi-Shaker  Test  and  Control  IEST-RP- 
DTE022.18  and  Welch’s  method.  In  particular,  an  adequate  number  of  statistical  degrees  of 
freedom  (DOFs)  have  been  obtained  to  provide  infonnation  with  acceptable  statistical  error. 
Consideration  must  be  given  to  not  only  statistical  error  in  auto-spectral  density  estimates,  but 
also  in  cross-spectral  density  estimates  (including  transfer  and  coherence  function  estimates). 

4.3.1  Instrumentation. 

For  the  purpose  of  this  TOP,  all  instrumentation  related  discussions  will  be  limited  to  linear 
accelerometers  and  engineering  units  of  g’s,  as  was  the  case  in  the  general  control  discussion 
provided  in  paragraph  3.1.1.  Linear  accelerometers  have  several  advantages  including 
familiarity  to  most  users,  low  cost,  wide  bandwidth,  small  size  and  weight,  and  readily  available 
low  cost  highly  reliable  signal  conditioning  options. 
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4.4  Use  of  Rigid  Body  Modes. 

a.  In  single  axis  testing,  the  control  input  is  often  defined  with  a  single  accelerometer. 
This  is  satisfactory  if  the  shaker  and  test  fixtures  are  rigid  within  the  frequency  band  of  interest. 

If  the  shaker  and  test  fixtures  are  not  rigid,  the  technique  of  using  a  single  accelerometer  for 
control  can  sometimes  lead  to  serious  difficulty.  To  overcome  these  problems,  methods  using 
the  average  of  several  accelerometers  and/or  force  limiting  have  come  into  common  practice.  In 
MEMA  testing,  the  problem  can  be  more  serious  as  non-rigid  body  response  is  more  common. 
When  considering  the  special  case  of  multiple  shakers  exciting  a  test  item  with  multiple  rigid 
body  degrees  of  freedom,  the  use  of  the  input  transfonnation  to  define  the  response  in  terms  of 
rigid  body  modes  has  several  advantages.  It  is  somewhat  analogous  to  a  generalization  of  the 
common  practice  for  single  axis  testing.  If  there  are  more  control  channels  than  rigid  body 
degrees  of  freedom,  and  an  input  transfonnation  matrix  is  defined  to  transfonn  the  control 
accelerometers  into  rigid  body  modes,  one  essentially  defines  the  motion  of  each  rigid  body 
mode  as  a  weighted  average  of  the  accelerometers  active  for  the  mode.  In  many  cases,  given  the 
control  authority  of  the  shakers,  this  is  about  the  best  viable  solution.  It  is  analogous  to 
averaging  accelerometers  for  a  single  axis  test,  which  is  common  practice.  The  elastic  modes  are 
not  controlled,  since  often  the  control  authority  over  these  modes  does  not  exist.  The  system  is 
driven  with  an  equivalent  rigid  body  motion  in  each  of  the  rigid  body  modes.  It  is  necessary  to 
make  sure  that  for  any  mode  the  transformation  of  the  control  accelerometers  {a}  does  not  result 
in  zero  for  any  of  the  rigid  body  modes.  If  higher  flexural  modes  are  present  they  will  not  be 
controlled.  In  theory  the  flexural  modes  can  be  controlled  by  adding  control  variables,  but  this 
requires  knowledge  of  the  modes  in  the  test  setup.  This  infonnation  can  only  be  determined  with 
materiel  in  the  test  configuration.  For  this  reason,  it  is  sometimes  desirable  to  allow  modification 
of  the  test  requirements  after  this  information  is  made  available.  Exactly  how  this  will  be 
accomplished  in  specification  writing  will  have  to  be  detennined  at  a  later  date. 

b.  An  advantage  of  using  rigid  body  modes  in  the  specification  is  that  the  field 
measurements  used  to  define  the  environment  can  be  made  with  the  transducers  in  locations 
different  from  the  locations  of  the  transducers  used  in  the  laboratory  test.  The  field 
measurements  are  reduced  to  equivalent  rigid  body  modes  using  an  acceleration  transfonnation 
matrix  (refer  to  paragraph  4.4.1),  and  the  modes  are  controlled  on  the  test  using  another 
transfonnation  matrix  for  the  laboratory  test  configuration.  The  two  transfonnation  matrices  do 
not  have  to  be  the  same.  Use  of  alternate  control  points,  while  maintaining  a  full  rank 
transfonnation  matrix,  provides  a  way  of  making  the  laboratory  test  “equivalent”  in  the  sense  of 
the  rigid  body  modes. 

c.  A  practical  difficulty  arises  when  more  modes  are  attempted  to  be  controlled.  The 
general  case  of  six  (6)  rigid  body  modes  requires  the  specification  of  a  6  x  6  SDM  (6  ASD’s  and 
15  USD’s).  Physical  understanding  of  the  SDM  matrix  associated  with  rigid-body  motion  by 
itself  is  difficult  without  the  additional  complications  of  elastic  DOFs.  Furthennore,  it  is  difficult 
to  assure  that  the  specification  results  in  a  positive  definite  SDM,  which  is  a  physical 
requirement.  (Additional  discussion  on  positive  definite  matrices  is  the  subject  of  paragraph 
4.5.2.) 
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4.4. 1  Acceleration  (Input)  Transfonnation. 


The  acceleration  to  control  space  transfonnation  matrix,  Tu  ,  commonly  referred  to  as  the  “input 
transformation  matrix”  from  the  control  system  perspective,  is  defined  in  the  article  “Applying 
Coordinate  Trans fonnations  to  Multi-DOF  Shaker  Control”9  and  generalized  in  the  article 
“Benefits  and  Challenges  of  Over-Actuated  Excitation  Systems”10.  The  acceleration 
transformation  matrix  transforms  a  set  of  accelerometer  measurements  into  a  set  of  control 
variables.  Often  these  control  variables  are  descriptions  of  rigid  body  modes.  The  acceleration 
transformation  is  usually  performed  in  the  time  domain  as: 


M=t,W 


4.4. 1 . 1  Acceleration  (Input)  Transformation  Derivation. 


One  goal  of  this  TOP  is  to  define  a  standard  nomenclature.  The  following  summary  has  been 
restructured  to  the  nomenclature  defined  by  this  TOP.  Referring  to  the  input  transformation 

derivation10,  a  generic  acceleration  measurement  at  the  k,h  position  in  orientation  j  is  structured 
as  Equation  4.1: 


(4.1) 


where  a0  is  the  linear  acceleration  at  some  reference  point  designated  the  “origin”,  a  is  the 
angular  acceleration  of  the  body  (assuming  it  is  rigid),  k  e  (1, 2 . Na ) ,  /  e  (l,  2 . n*) ,  j  e  (x,y,z) , 

and  =  [1  0  0] ,  e^y  =  [0  10],  and  eTz  =  [0  0  l]  are  row  selection  vectors  (as  shown  assuming 
accelerometer  orientation  is  aligned  per  a  traditional  right  hand  Cartesian  system).  Parameter  Na 

represents  the  number  of  accelerometer  measurements  (as  previously  defined)  and  n*  <Na  the 

number  of  measurement  locations;  e.g.,  utilization  of  multi-axis  accelerometers  results  in 

n*  <  Na  .  Vector  r  is  the  position  vector  relating  the  position  of  measurement  location  i  to  a  user 

”  p  p~  x 

defined  origin.  r,  is  the  skew  symmetric  operator  equivalent  of  the  cross  product,  making 
the  matrix  based  computations  in  Equation  4. 1  possible.  The  matrix  equivalent  of  a  vector  (i.e., 
a  coordinatized  vector  quantity)  is  denoted  as  (  1  (_)'  j  where  the  right  superscript  and  subscript 
identify  the  body  and  point  of  interest  respectively,  and  the  left  superscript  denotes  the 

p  p 

coordinate  frame  in  which  the  vector  quantity  was  coordinatized;  e.g.,  r,  in  Equation  4. 1 
denotes  the  ith  point  on  body  P  (the  platform)  coordinatized  in  frame  J~P  -  the  platfonn’s 
coordinate  frame. 
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4.4. 1 .2  Equation  4. 1  represents  one  equation  in  six  unknowns,  the  three  components  of  the  linear 
acceleration  of  the  reference  point  and  the  three  components  of  the  rigid  body  angular 
acceleration.  In  order  to  determine  these  quantities,  at  least  six  measurements  are  needed.  These 
requirements  are  not  as  stringent  as  that  reported  in  the  article  “On  the  Use  of  Linear 
Accelerometers  in  Six -DOF  Laboratory  Motion  Replication”11  because  of  the  assumptions  above 
(i.e.,  small  angular  velocities  and  rigid  body). 

Let’s  consider  the  most  general  case  of  Na  measurements  from  n*  locations.  In  this  case, 
Equation  4.1  becomes: 


T  T 

rp  pi 

X 

\ " 

e-j  ~e-j 

* — t 

r 

a2, 

e]  -e] 

rp  p i 
r; 

X 

mm 

akj 

— 

J 

— 

-J  -J 

L  >  J 

p  p 

a 

&ni  _ 

T  T 

6  /  -e  : 

_p  p 

r\ 

X 

(6x1) 

(nxl) 

— J  -J 

n 

/e(l,2 . /)*),/'  e(x,y,z) 


(nx6) 


which  using  the  nomenclature  defined  in  this  TOP  is  of  the  form: 


Mm 

l  >  Meas 
(nxl) 


Motion 


(nx6)  (6x1) 


(4.2) 


where  {c}Motion  is  a  6  x  1  matrix  of  unknown  linear  and  angular  accelerations  and  {a}  is  an 
nxl  matrix  of  acceleration  measurements.  Observe  that  [Va  J  is  entirely  defined  by  knowledge  of 
(i)  placement,  (ii)  orientation,  and  (iii)  utilized  signals  of  the  accelerometers. 


Observe  that  if  Ta  is  of  full  column  rank,  then 
follows: 


T  T 

a  a 


exists  enabling  {c}  to  be  solved  as 


Jai  =]’  c 

t  /  Meas  a  \  /  Motion 

t'  M„«„  =  T,TT,  i  c  !•  . . . 


[TXT  T;(a}Mm=[T;T,J  t;t.{c}Moi1o 
[T.TT,J'T,qa}Mm={c}, 


'  Motion 


Defining  Ta  =  TaTTa  TaT ,  Equation  4.2  can  be  rewritten  as: 


Motion 


Meas 


(4.3) 
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Where  [T,]  is  a  6  x  n  matrix  referred  to  in  the  literature  as  the  “Acceleration  Transfonn  Matrix’ 


TtT 

a  a 


exists  in  order 


or  “Input  Transform  Matrix”.  Observe  that  the  critical  requirement  that 

to  derive  the  input  transformation  matrix  [Ta] ,  is  solely  a  function  of  placement  and  orientation 
of  measurement  transducers. 


4.4.2  Drive  (Output)  Transformation. 

a.  Although  details  of  the  Drive  Transformation  are  not  required  to  develop  a  MDOF 
VSD  reference,  a  short  summary  of  the  concept  is  provided  for  general  knowledge.  Referring  to 
the  schematic  in  Figure  2,  transfonnation  matrix  Ts  transforms  the  Nd  drive  variables  into  Ns 

shaker  drive  signals.  Reference  10  provides  a  formal  derivation  of  the  transformation  matrix,  Ts. 
Note  that  while  the  “acceleration  transfonnation”  was  computed  based  on  knowledge  of  position 
and  polarity  of  the  control  accelerometers,  the  transformation  matrix,  Ts  is  dependent  upon  the 
position  and  line  of  action  (LOA)  of  the  individual  actuators.  In  this  TOP  and  within 
reference  10  Ts  is  refened  to  as  the  “drive  transformation”  or  “output  transfonnation”.  The 
following  cases  summarize  the  computation  of  Ts  and  the  effect  on  the  control  process. 

(1)  Case  1 :  Configurations  in  which  the  number  of  motion  degrees-of-freedom  or 
control  signals,  Nc  and  the  number  of  output  control  variables,  Nd  are  the  same  is  refened  to  as 
“square”  control.  If  the  number  of  output  control  variables,  Nd  and  the  number  of  shakers,  Ns  is 
the  same,  the  transformation  matrix,  Ts  will  simply  be  the  Identity  matrix. 

(2)  Case  2:  Configurations  in  which  the  number  of  shakers  Ns  exceeds  the  number  of 
output  control  variables  Nd  ,  the  excitation  system  is  said  to  be  over-detennined  or  over- actuated. 
In  such  cases,  some  of  the  drives  will  be  linear  combinations  of  other  drives.  Furthennore,  if  Ts 
is  a  constant  which  is  employed  in  the  time  domain,  the  individual  actuators  must  be  matched 
(e.g.  matched  frequency  response  functions  (FRFs)). 

(3)  Case  3:  Configurations  in  which  the  number  of  shakers,  Ns  is  less  than  the 
number  of  control  signals,  Nc ,  the  excitation  system  is  said  to  be  under-determined  or  under¬ 
actuated.  In  such  cases,  exact  control  of  the  SDM  is  not  possible. 

b.  In  theory,  Ts  could  be  implemented  before  or  after  the  transformation  into  the  time 
domain.  One  advantage  of  placing  the  transformation  in  the  frequency  domain  section  of  the 
control  algorithm  is  that  the  matrix  could  then  be  made  a  function  of  frequency.  Having  the 
transformation  matrix,  Ts,  a  constant  assumes  the  shakers  are  matched  and  the  desired 
transformation  can  be  deduced. 
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4.4.2. 1  Drive  (Output)  Transformation  Derivation. 

a.  As  previously  stated,  one  goal  of  this  TOP  is  to  recommend  a  standard  nomenclature. 
The  following  summary  from  reference  number  10  has  been  restructured  to  the  nomenclature 
recommended  by  this  TOP.  Figure  3  illustrates  the  generalized  multi-axis  vibration  system. 


Figure  3.  Generalized  multi-axis  vibration  system. 


b.  Refer  to  reference  number  10  for  a  detailed  derivation  of  Equation  4.4.  The  following 
summary  illustrates  how  the  output  transform,  Ts  is  associated  with  the  P-Matrix,  (Plucker 
Matrix)  discussed  in  the  reference. 
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F 

N,x  1 


c 
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(4.4) 


c.  In  Equation  4.4,  [P]  represents  the  Plucker  Matrix  which  is  derived  from  known 
geometric  parameters  associated  with  the  individual  actuators,  [F]  represents  the  drive  and  [C] 
represents  the  desired  motion.  The  variables  11  lc  represent  the  LOA  vectors  for  each  of  the 
actuators  and  Pmf  is  the  moment  arm  associated  with  force  f] .  Observe  that  the  maximum 
dimension  for  the  [C]  matrix  will  be  six,  if  all  six  traditional  motion  DOFs  are  being  considered 
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(i.e.  Nd  =6).  As  stated  in  paragraph  4.4.2,  Case  1  scenarios  will  simply  have  an  identity  matrix 
as  the  output  transformation  matrix  and  Case  3  scenarios  (under-actuated)  will  not  have  a  unique 
solution.  Case  2  scenarios  (over-actuated)  may  be  addressed  in  terms  of  output  transformations. 
The  objective  is  to  determine  [F]  in  Equation  4.4,  yielding  the  Ns  drive  signals  as  follows: 

(1)  Define  F  =  P  1  D  and  substitute  into  P  F  =  C  yielding  P  P  1  D  =  C 

Nsx  1  Nsx6  6x1  6 xNs  Nsx  1  6x1  6xJVs  Nsx 6  6x1  6x1 

(2)  P  P  will  be  of  full  rank  (i.e.  invertible)  if  [ Pi  is  of  full  rank. 

6  xNs  Nsx  6  L  J 


(a)  If  [P]  is  of  full  rank: 


D 

6x1 


"1-1 


P  P  ' 

6  xN,  N.x  6 


c 

6x1 


(b)  If  [P]  is  not  full  rank,  actuator  placement  is  not  sufficient  to  obtain  the 
mechanical  DOF’s  desired. 


(3) 


Substituting  results  from  (2)  yields 


F  =  P  r  D  =  Pr  P  Pr 

Nsxl  Nsx6  6x1  Nsx6  6xNs  Nsx6 


C 

6x1 


(4) 


T  =  P  7  P  P  1 

Nsx  6  Nsx6  L  6xNs  Nsx6 


d.  The  discussions  within  this  paragraph  and  previous  derivation  assumed  Nd  -  6  .  In  the 
event Nd  <6,  Nd  would  represent  the  actual  number  of  mechanical  DOFs.  In  terms  of  the 
nomenclature  of  Figure  2,  and  assuming  matched  actuators  are  employed,  voltage  drives  to  the 
shakers  for  the  over- actuated  scenario  would  be  defined  as  Equation  4.5: 

{s(  =  T,{d}  (4.5) 


4.5  Data  Analysis. 

a.  Ensure  transducer  placements  have  been  addressed,  to  guarantee  the  desired  motion 
DOFs  may  be  resolved  (refer  to  paragraph  4.4. 1.2),  and  that  common  data  validity  checks  are 
performed.  Then,  it  is  recommended  that  appropriate  combinations  of  the  linear  acceleration 
measurements  be  transformed  into  the  6  traditional  motion  DOFs  through  implementation  of  the 
acceleration  transformation  matrix.  The  transformed  time  histories  will  be  referenced  to  a  single 
point  on  the  structure  referred  to  as  the  “origin”  as  discussed  in  paragraph  4.4.1. 

b.  A  SDM  for  each  test  configuration  identified  in  the  mission  scenario  should  be 
computed.  In  addressing  the  VSD  techniques  for  reducing  an  ensemble  of  data,  in  this  case  an 
ensemble  of  SDM’s,  the  analyst  will  be  required  to  deal  with  the  ASD  terms  (the  diagonal  terms 
of  the  SDM)  and  CSD  tenns  (the  off-diagonal  terms  of  the  SDM). 
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4.5.1  Phase  and  Coherence  based  Representations  of  CSD  Terms. 


Although  the  off-diagonal  terms  of  the  SDM  are  computed  in  tenns  of  a  CSD,  it  is  common 
among  control  system  vendors  to  allow  cross  tenns  to  be  defined  in  terms  of  Phase  and 
Coherence.  This  is  a  convenient  option  in  that  it  is  often  easier  to  physically  interpret  SDM  CSD 
terms  in  terms  of  Phase  and  Coherence.  There  is  a  direct  relationship  between  the  two 
techniques  of  defining  the  cross  tenns  of  the  SDM  that  is  dependent  upon  the  definition  of 

I  |2 


ordinary  coherence  between  two  signals, 


y2  = 


G„G, 


Normalizing  the  CSD  tenns  of  the  SDM 


by  yields  a  normalized  spectral  density  matrix  (SDMn)  in  which  the  ASD  terms  are  not 


affected  and  the  magnitude  of  the  normalized  CSD  terms  are  defined  as 


,  which  is 


equivalent  to  the  square  root  of  the  ordinary  coherence  function,  while  not  affecting  the  original 
phase  relationship  of  the  CSD  terms.  Similarly,  the  normalized  spectral  density  matrix,  SDMn, 
may  be  transformed  back  to  the  original  CSD  form  of  the  SDM. 


4.5.2  Positive  Definite  SDM  Considerations. 


a.  Any  specified  spectral  density  matrix  must  be  positive  semi-definite  to  be  physically 
realizable.  In  practice  it  must  be  positive  definite.  The  detenninate  of  the  matrix  must  be  >0. 
All  the  eigenvalues  of  the  SDM  must  be  >0.  This  must  be  true  at  all  frequencies.  It  must  be 
possible  to  perform  a  Cholesky  decomposition  of  the  specified  SDM.  Another  property  of 

j2 

positive  semi  definite  matrices  is  from  Matrix  Computations  : 

i  i2 

2  ®J 

® ,  <  ® ..®  or  0  <  y2  =  <  1 

!/  ••  jj  '  00 

«  JJ 


In  the  terms  of  random  vibrations  the  ordinary  coherence,  y1  between  signals  must  be  less  than 
or  equal  to  one.  In  practical  terms,  if  the  coherence  between  any  pair  of  signals  is  one,  the  SDM 
will  be  positive  semi-definite  and  the  control  system  will  have  problems.  Note  that  in  general,  if 
D  is  Hennitian  and  positive  semi-definite  C  will  also  be  Hermitian  and  positive  semi-definite. 

b.  If  all  the  eigenvalues  are  non-negative,  the  matrix  is  positive  semi-definite.  If  any  of 
the  eigenvalues  are  zero,  it  implies  that  one  or  more  of  the  rows  of  the  spectral  density  matrix  are 
a  linear  combination  of  other  rows.  In  practice,  one  would  typically  expect  to  deal  only  with 
positive  definite  matrices.  Observe  that  even  a  small  amount  of  noise  or  nonlinearity  will  result 
in  a  positive  definite  matrix.  If  a  matrix  is  positive  definite,  the  matrix  can  always  be  factored 
using  Cholesky  decomposition, 


<D  =  LL' 


where  L  is  a  lower  triangular  matrix.  Which  without  loss  of  generality  can  be  rewritten  as, 
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O  =  LIL' 

where  I  is  the  identity  matrix.  In  this  application,  I  is  not  really  the  identity  matrix.  I  is  a 
spectral  density  matrix.  At  every  frequency,  I  is  a  diagonal  matrix  of  ones.  The  components  in  I 
are  independent  since  all  the  off  diagonal  elements  are  zero.  It  is  now  clear  why  the  cross 
spectral  density  matrix  must  be  positive  definite.  If  any  of  the  elements  in  I  are  zero,  it  implies 
that  there  are  less  than  N  (the  number  of  rows  or  columns  in  0 )  independent  sources  in  0 . 

Some  of  the  rows  and  columns  are  linear  combinations  of  other  rows  and  columns.  The  identity 
matrix  is  positive  definite,  therefore  0  must  be  positive  definite.  Using  the  interpretation  of 
Random  Data  Analysis  and  Measurement  Procedures13,  the  diagonal  elements  of  I  can  be 
interpreted  as  the  auto-spectral  densities  of  independent  random  noise  sources.  The  maximum 
number  of  independent  noise  sources  is  N.  If  some  of  the  elements  in  I  are  zero,  the  problem  can 
still  be  solved  by  making  the  corresponding  rows  and  columns  of  L  zero.  This  is  the  positive 
semi-definite  case.  This  case  corresponds  to  the  case  where  there  exists  less  than  N  independent 
sources.  Some  of  the  N  sources  are  linear  combinations  of  other  sources.  This  case  will  be  very 
difficult  for  the  control  system.  In  general  one  may  make  some  of  the  sources  small  but  not  zero. 
Part  of  this  document  will  discuss  the  generation  of  a  desired  control  SDM  to  make  the  control 
problem  achievable  and  hopefully  relatively  easy  for  the  control  system  to  implement. 

c.  In  general  the  control  problem  is  an  inverse  problem.  The  desired  control  SDM  (the 
output  of  the  system  under  test)  is  known,  and  the  drive  (input  to  the  system  under  test)  SDM 
must  be  computed.  There  is  a  potential  point  of  confusion  here.  The  control  system 
manufacturers  think  the  drive  SDM  is  the  output  of  the  control  system,  which  is  the  input  to  the 
shaker  system.  Similarly,  the  control  system  input  is  the  output  of  the  shaker  system. 

Paragraphs  4.4.1  and  4.4.2  provide  nomenclature  employed  for  input  and  output  transformations 
as  they  are  applied  within  this  document. 

d.  Inverse  problems  can  be  very  difficult  as  multiplication  by  a  matrix  inverse  is 
required.  If  the  matrix  is  ill-conditioned,  the  result  will  be  similar  to  dividing  by  zero  for  the 
scalar  case. 

For  the  case  in  which  the  number  of  inputs  and  outputs  are  the  same;  H  is  a  square  matrix  of 
FRF’s.  The  solution  is  to  invert  H.  The  solution  for  the  drive  matrix  is  then  given  by: 

Z  =  H  1 
D  =  ZRZ' 

This  of  course  assumes  H  is  well  conditioned  and  the  inverse  exists.  Part  of  this  document  will 
discuss  issues  to  help  the  process  of  achieving  a  well  conditioned  H  matrix. 

The  H  matrix  is  typically  estimated  from: 


H  =  ScdD 
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The  inverse  of  D  must  exist.  This  implies  that  D  must  be  positive  definite.  The  initial  estimate 
of  H  is  determined  by  exciting  the  system  with  a  set  of  independent  white  inputs  in  a  pretest 
environment.  If  H  is  to  be  corrected  during  the  test,  D  must  be  positive  definite  during  the  test 

or  special  provisions  must  be  used  to  avoid  the  inversion  of  D  at  frequencies  where  D  is  not 
positive  definite.  This  is  one  of  the  reasons  the  reference  R  rarely  has  any  of  the  coherences 
equal  to  unity. 

4.5.3  Data  Compression. 

a.  Use  of  time  compression  techniques  such  as  Miner-Palmgren  may  be  employed  to 
modify  the  ASD  terms.  References  numbers  1  and  3  provide  discussions  on  time  compression. 

In  the  simplest  terms,  the  Miner-Palmgren  Hypothesis  (Miner’s  rule)  is  a  set  of  mathematical 
equations  used  to  scale  vibration  spectra  levels  and  their  associated  test  times.  It  provides  a 
convenient  means  to  analyze  fatigue  damage  resulting  from  cyclical  stressing.  The  mathematical 
expression  and  variable  descriptions  for  this  technique  are  illustrated  below  in  Equation  4.6: 

t_2_  = 

h 

where: 

6  =  equivalent  test  time 

t2  =  in-service  time  for  specified  condition 

=  severity  (root  mean  square  ((rms))  at  test  condition 

S2  =  severity  (nns)  at  in-service  condition 

(The  ratio  SjS2  is  commonly  known  as  the  exaggeration  factor.) 

M  =  a  value  based  on  (but  not  equal  to)  the  slope  of  the  S-N  curve  for  the  appropriate 
material  where  S  represents  the  stress  amplitude  and  N  represents  the  mean  number  of 
constant  amplitude  load  applications  expected  to  cause  failure.  For  the  MDOF  VSD  work 
at  hand,  the  default  of  M  =  7  was  selected  per  reference  number  1 . 

b.  It  is  recommended  that  the  final  vibration  specification  ASD  terms  are  no  greater  than 
3  decibel  (dB)  higher  than  maximum  values  measured  in  the  field.  Miner-Palmgren  will  be 
employed  to  the  ASD  portion  of  the  SDM  in  the  same  manner  as  one  would  employ  for  a 
traditional  1-DOF  scenario.  Details  such  as  maintain  common  test  durations  between 
mechanical  DOFs  are  address  in  Paragraph  6. 

4.5.4  Limiting  Strategies. 

Traditional  notching  techniques  may  also  be  employed  if  impedance  mismatches  lead  to 
unrealistically  high  test  item  response.  Notching  techniques  may  be  employed  across  all 
actuators  with  equal  weighting  or  by  weighting  notching  at  each  actuator  as  a  function  of 
coherence  between  the  actuators  and  the  location  of  interest.  In  addition  to  traditional  notching 
based  on  acceleration  spectra,  it  is  also  possible  to  consider  limiting  based  on  other  parameters 


(4.6) 
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(e.g.  von  Mises  Stress  or  Force  limiting).  As  with  any  notching  scheme,  it  is  critical  that  any 
resulting  deviations  to  the  test  or  test  tolerances  must  be  approved  by  the  appropriate  test 
authority  and  must  be  clearly  documented  in  the  test  plan  and  final  report. 

4.5.5  Minimum  Drive  Considerations. 

A  number  of  challenges  have  been  identified  in  addressing  the  objective  of  establishing  a 
reference  SDM  for  multiple  exciter  test  (MET)  scenarios.  One  major  area  of  concern  is  related 
to  the  fact  that  it  is  highly  likely  that  there  will  be  mechanical  impedance  differences  between  the 
field  and  laboratory  conditions.  Given  these  impendence  mismatch  issues,  it  is  undesirable  to 
force  the  test  item  into  what  could  potentially  be  an  unnatural  state  as  fixtured  in  the  laboratory. 
Optimally,  achieving  the  specified  autospectra  without  excessively  taxing  the  excitation  system 
is  desired.  Smallwood  made  a  general  approach  to  establishing  minimum  drive  criteria  in  the 
article  “MIMO  Linear  Systems  Extreme  Inputs/Outputs”14.  Unfortunately,  the  technique  does 
not  always  guarantee  the  resulting  SDM  to  be  positive  semi-definite. 

4.5.5. 1  Independent  Drives. 

a.  Although  an  active  area  of  research,  general  techniques  to  address  minimum  drive 
criteria  have  not  been  formally  established  at  the  time  of  this  publication.  A  proposed  approach 
for  trending  drive  voltages  towards  minimums  while  maintaining  a  positive-  definite  SDM,  is 
discussed  in  the  article  “A  Proposed  Method  to  Generate  a  Spectral  Density  Matrix  for  a  MIMO 
Vibration  Test”15,  and  is  summarized  below: 

(1)  Taking  a  clue  from  the  modal  test  community,  assume  the  drive  signals  to  the 
excitation  system  will  be  uncorrelated.  Typically  for  a  vibration  test,  the  drives  are  the  voltage 
inputs  to  the  shakers.  For  a  simulation,  the  inputs  into  a  model  are  often  forces.  It  is  always 
possible  to  excite  the  system  with  uncorrelated  inputs.  This  is  standard  practice  in  the  modal 
community,  and  is  standard  practice  when  perfonning  the  system  identification  for  MIMO  test 
systems.  This  leads  to  the  logical  question:  Is  it  possible  to  generate  a  set  of  uncorrelated  inputs 
that  will  produce  a  desired  set  of  response  autospectra  (the  diagonal  of  the  output  SDM)? 

(2)  The  general  equation  relating  the  control  point  accelerations  to  the  drive  voltages 
is  given  in  Random  Vibrations,  Theory  and  Practice16: 

SY  =  HSXH' 

where  H'  is  the  conjugate  transpose  of  H ,  and  Sx  and  Sy  are  SDM’s.  H  is  a  matrix  of 
frequency  response  functions  relating  the  output  to  the  input  of  the  excitation  system.  In  our 
case,  ideally,  Sx  will  be  a  diagonal  matrix.  Let  X  be  a  column  vector  of  the  diagonal  of  Sx  or, 

X  =  diag(Sx) ,  and  Y  =  diag(SY) .  The  relationship  between  the  autospectra,  as  shown  in 
Annex  D  proof  1 ,  is  given  by: 


Y  =  HX 
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where: 


H  =  H.*conj(H) 

where:  *  indicates  an  element  by  element  multiplication.  Htj  =|  Hfj  |2 . 

The  solution  is  given  by: 

\  TT  v 

b.  In  some  cases  the  result  will  include  negative  elements  inX  .  This  is  not  physically 
possible.  It  indicates  that  the  desired  ASD’s  cannot  be  achieved  with  independent  drives.  In  this 
case  the  negative  values  are  set  to  zero,  and  the  output  SDM  is  recomputed  from  SY  =  HSXH' 
using  the  modified  input  spectral  density  matrix  (the  negative  values  set  to  zero).  The  resulting 
control  point  acceleration  autospectra,  will  not  be  at  the  desired  levels.  To  correct  this  problem, 
the  control  point  acceleration  autospectra  are  rescaled  to  the  desired  levels,  keeping  the  phase 
and  ordinary  coherence  the  same.  This  is  accomplished  by  pre  and  post  multiplying  the  SDM  by 
a  diagonal  matrix  whose  elements  are  the  square  root  of  the  ratio  of  the  desired  ASD  to  the 
computed  ASD: 


C  _  c  c  c 

^Ynew  ^  v 15  Void  °  s 


where  Ss  is  a  diagonal  matrix  and 


S. ,  = 


l, old 


Note:  Setting  S Yii,new=  1 ,  provides  an  efficient  way  to  compute  the  normalized  SDM  where  the 
diagonals  are  one  and  the  magnitude  of  the  off  diagonals  squared  are  the  ordinary  coherence  and 
the  phase  of  the  off  diagonal  elements  is  the  phase  of  the  cross  spectra. 

The  drive  SDM  can  then  be  computed  as: 


Sxnew  -  ZSynewZ’ 


where  Z  =  pinv(H) ,  the  Moore-Penrose  pseudo  inverse.  If  H  is  square  and  full  rank,  the 
solution  typically  ends  here.  If  H  is  not  square  or  not  full  ranked: 

s¥™=  HS^H' 

Thediag(SYm;H)  may  not  yield  the  desired  ASD’s.  In  this  case,  an  iterative  approach  will  often 
improve  the  result. 
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4.6  Independent  References. 

a.  It  is  sometimes  desirable  to  define  the  reference  spectrum  in  terms  of  a  diagonal 
matrix  of  autospectra.  Several  reasons  drive  us  in  this  direction.  One  common  case  is  that  only 
knowledge  of  the  autospectra  from  the  field  environments  is  available.  Several  factors  can  result 
in  this  situation.  First  the  field  data  may  have  been  acquired  without  phase  information.  Second, 
the  resulting  cross  spectra  can  have  a  very  complicated  structure  which  is  impractical  to 
implement  in  a  specification.  Enveloping  amplitudes  is  possible,  but  enveloping  the  phase  is 
much  more  difficult.  Third,  the  specification  may  be  a  composite  of  several  environments, 
making  the  definition  of  cross  spectra  very  difficult.  Fourth,  the  vehicle  on  which  the  field  data 
were  taken  may  not  be  identical  to  the  test  vehicle.  Fifth,  the  boundary  conditions  in  the  field 
may  be  different  from  the  boundary  conditions  in  the  laboratory. 

b.  Small  changes  in  the  modal  frequencies  caused  by  any  of  the  above  factors  can  change 
the  phase  at  any  frequency  near  a  mode  by  a  large  amount.  All  these  factors  make  the 
specification  of  the  cross  spectra  difficult.  An  option  is  to  ignore  the  cross  spectra  and  set  them 
all  to  zero.  This  has  the  theoretical  advantage  of  providing  an  excitation  that  in  some  sense 
covers  the  control  variable  response  space. 

c.  The  drive  signals  can  readily  be  computed  yielding  uncorrelated  motion  (in  this  case 
the  SDM  of  the  uncorrelated  reference  spectra  Y  is  diagonal)  from: 

SX0=ZSY0Z’ 

This  approach  is  currently  available  in  commercial  control  systems.  You  simply  specify  the 
reference  SDM  as  a  diagonal  matrix  with  the  cross  spectra  (or  equivalently  the  coherences)  zero 
or  near  zero.  This  is  typically  a  conservative  approach. 

d.  In  contrast  to  the  independent  drive  discussion  in  paragraph  4.5.5. 1,  the  danger  with 
the  independent  reference  concept  is  that  this  specification  of  control  variables  may  be  overly 
conservative  near  frequencies  dominated  by  a  single  mode.  An  important  clue  that  the  result 
may  be  overly  conservative  is  the  trace  of  the  drive  voltages.  This  trace  should  be  monitored  and 
if  overly  large  in  some  band  of  frequencies,  limits  can  be  negotiated  and  implemented. 

4.7  Recommended  Practices  Summary. 

The  following  list  provides  recommendations  and  general  guidance  to  be  considered  when 
addressing  the  multi-axis  VSD. 

a.  If  possible,  specify  the  test  in  terms  of  the  rigid  body  motion. 

b.  Over  specification  of  the  control  accelerometers  is  desirable.  Use  more  control 
accelerometers  than  degrees  of  freedom  in  the  test. 
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c.  If  possible,  the  entire  SDM  should  be  specified.  A  method  to  automate  the  generation 
of  envelopes  may  be  desired.  This  will  permit  the  generation  of  the  envelopes  to  be  less 
developer  specific. 

d.  If  the  entire  SDM  is  specified,  it  is  suggested  that  the  coherence  be  set  to  near  zero  if 
the  desired  coherence  is  below  0.2.  It  should  be  recognized  that  the  estimation  of  coherence  is  a 
biased  result  (the  result  will  always  be  positive).  It  is  recognized  that  the  estimated  coherence 
will  never  be  zero;  however,  the  control  software  can  attempt  to  make  the  coherence  as  low  as 
possible.  The  tolerance  on  the  coherence  must  recognize  the  bias.  If  the  coherence  is  small  the 
phase  is  not  important.  For  convenience,  establishing  a  zero  phase  is  a  reasonable  specification 
when  the  coherence  is  low. 

e.  If  step  c  becomes  too  complicated,  it  is  recommended  that  the  test  be  run  with  near 
zero  coherence. 

f.  If  step  e  results  in  unrealistic  responses,  try  using  the  independent  drive  option. 

g.  Consider  a  compromise  position  between  independent  reference  criteria  of  step  e  and 
independent  drive  criteria  as  recommended  in  step  f. 

h.  If  the  drive  requirements  are  excessive  at  some  frequencies,  allow  the  test  to  be 
modified  to  reduce  the  drive  requirements  as  discussed  in  paragraphs  4.5.4  and  4.5.5. 

i.  It  is  understood  that  MIMO  testing  is  more  complicated  than  single-input  single  output 
(SISO)  testing.  The  specifications  must  reflect  the  desires  of  a  knowledgeable  environmental 
test  engineer.  Good  communication  between  the  project  team,  the  environmental  test  engineer 
and  the  test  lab  must  be  maintained  to  achieve  the  desired  test  results. 

5.  DATA  REQUIRED. 

•5 

Just  as  discussed  in  NATO  AECTP  Leaflet  2410  ,  field  data  must  be  acquired  based  upon  the 
anticipated  mission  scenario  of  the  unit  under  test  (UUT).  As  detailed  in  paragraph  4.4. 1 . 1  and 
reference  number  1,  transducer  placement  and  orientation  are  critical  and  must  be  thoroughly 
documented. 

5.1  Reference  SDM  Development. 

As  stated  in  paragraph  4.5,  a  SDM  in  terms  of  the  six  traditional  rigid  body  modes  should  be 
computed  for  each  test  configuration  identified  in  the  mission  scenario. 

5.1.1  SDM  Ensemble  CSD  Characteristics. 

Based  on  the  characteristics  of  the  CSD  terms  of  the  ensemble  of  SDMs,  the  VSD  process  will 
yield  a  vibration  specification  consistent  with  one  of  the  three  cases  that  follow: 
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a.  Case  1.  Coherence  Terms  Approaching  Zero  (Independent  Motion  DOFs)  -  This  is 
the  easiest  situation  to  deal  with  in  that  each  motion  DOF  ASD  may  be  addressed  individually 
via  the  same  techniques  employed  in  1-DOF  VSD  as  discussed  in  reference  number  3.  When 
programming  the  vibration  control  system,  it  is  recommended  that  coherence  be  set  to  a  low  non¬ 
zero  level  (i.e.  y2  =0.1)  over  the  test  bandwidth  of  interest.  For  such  a  small  coherence,  the 
phase  parameter  is  essentially  a  random  variable  and  establishing  a  phase  specification  is  not 
required. 

A  special  situation  that  may  lead  an  analyst  to  develop  a  MDOF  vibration  specification  with 
independent  motion  DOFs,  would  be  a  composite  specification  that  encompasses  multiple 
vehicles  (i.e.  a  composite  wheeled  vehicle  specifications  comparable  to  those  in  MIL-STD- 
810G,  Method  514).  As  each  vehicle  will  tend  to  have  its  own  CSD  characteristics,  it  is  not 
possible  to  define  CSD  terms  such  that  a  single  coherence  and  phase  relationship  addresses  each 
vehicle.  Enveloping  techniques  that  work  well  in  addressing  magnitude  based  ASD  terms  are 
simply  not  applicable  in  addressing  phase  relationships  between  mechanical  DOFs. 

b.  Case  2.  Non-Zero  Coherence  across  a  Portion  of  the  Test  Bandwidth  -  When 
developing  a  MDOF  vibration  specification  based  on  a  single  platfonn,  one  would  expect  the 
CSD  terms  measured  across  the  range  of  scenarios  addressed  in  the  mission  scenario  to  be 
similar  in  nature.  The  dynamic  characteristics  of  the  structure  and  often  the  proximity  of  the 
measurement  transducers  will  greatly  influence  the  CSD  characteristics.  There  are  often 
situations  in  which  coherence  between  motion  DOFs  are  high  and  phase  is  well  defined,  but  only 
over  a  portion  of  the  test  spectrum.  This  is  a  common  observation  on  many  wheeled  vehicles 
where  coherence  is  high  at  lower  frequencies  (i.e.  frequencies  below  50  Hertz  (Hz)  and  near  zero 
at  higher  frequencies.  In  such  scenarios,  one  would  only  establish  coherence  and  phase 
specifications  for  the  portion  of  the  spectrum  with  high  coherence.  The  remainder  of  the 
spectrum  would  be  treated  as  in  Case  1.  Also,  in  establishing  CSD  reference  criteria,  the  analyst 
must  ensure  the  resulting  criteria  is  physically  realizable  (refer  to  paragraph  4.5.3  for  additional 
detail). 


c.  Case  3.  Non-Zero  Coherence  across  the  Full  Test  Bandwidth  -  This  scenario  is 
comparable  to  Case  2  with  coherence  being  defined  across  the  entire  test  bandwidth.  It  is 
anticipated  that  this  would  be  the  least  likely  scenario  in  a  MDOF  VSD  effort.  However,  it  is 
also  the  configuration  that  will  be  the  most  difficult  to  deal  with  from  both  a  VSD  development 
aspect  and  from  an  implementation  perspective.  In  addition  to  the  issue  of  ensuring  the  resulting 
SDM  reference  is  physically  realizable,  the  classic  problem  of  mechanical  impedance  mismatch 
between  field  and  laboratory  are  often  major  concerns  in  implementing  a  fully  defined  SDM 
reference  criterion  for  a  laboratory  test.  Specifically,  if  the  mechanical  impedance  between  field 
and  laboratory  are  not  very  well  matched  (and  they  usually  are  not),  there  may  be  portions  of  the 
spectrum  in  which  coherence  may  be  significantly  different  than  specified  and  simply  not 
controllable.  While  this  situation  is  also  possible  in  Case  2,  it  is  almost  certain  to  be  an  issue  in  a 
scenario  such  as  Case  3,  in  which  the  entire  test  bandwidth  has  a  CSD  reference  criteria.  This 
topic  of  uncontrollable  coherence  associated  with  mechanical  impedance  mismatches  is  a  control 
issue  for  all  three  Cases  and  is  discussed  further  in  the  minimum  drive  consideration  of 
paragraph  4.5.6. 
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d.  Regardless  of  which  of  the  three  cases  the  SDM  is  characterized  by,  it  is  highly  likely 
that  there  will  be  mechanical  impedance  differences  between  the  field  and  laboratory  conditions. 
In  some  cases  these  impedance  differences  may  result  in  excessive  drive  signals.  Although  the 
various  control  system  vendors  address  this  situation  in  varying  degrees,  it  may  still  be  necessary 
to  address  this  issue  through  test  operator  intervention  via  techniques  such  as  those  identified  in 
paragraphs  4.5.4  and  4.5.5. 

5.2  Test  Tolerance  Recommendations. 

Setting  tolerances  for  a  MIMO  test  is  challenging  given  the  large  amount  of  infonnation 
encompassed  by  the  reference  autospectra  and  cross  spectra  involved.  Additionally,  the  overall 
energy  is  not  necessarily  distributed  evenly  about  each  mechanical  DOF  and  dominant  DOFs 
often  tend  to  dominate  the  control.  The  objective  here  is  to  establish  a  reasonable  starting  point 
in  establishing  test  tolerances.  Experience  with  specific  test  configurations  may  be  employed  to 
refine  the  basic  requirements  defined  below.  As  usual,  any  test  specific  test  tolerances  should  be 
clearly  documented  within  the  test  plan. 

a.  Autospectra(1):  ±3  dB  for  /  <500 Hz  and  ±6dB  for  /  >500 Hz. 

( 1  'The  portion  of  the  spectrum  that  actually  reaches  the  maximum  tolerance  limits  is 
anticipated  in  narrow  bandwidths.  The  tolerance  on  the  overall  Grms  level  of  each 
controlled  DOF  shall  be  within  +15%  of  the  corresponding  reference. 

b.  Cross  spectra:  Define  tolerances  in  terms  of  Phase  and  Coherence 

(1)  Coherence:  For  ordinary  coherence  in  the  range  0.5  <  y2  <  1.0 ,  set  the  tolerance 
to  be  ±0.1  (avoid  establishing  a  coherence  reference  or  tolerance  of  1.0). 

(2)  Phase:  If  y2  <  0.5  ,  any  phase  is  acceptable.  If  0.5  <  y2  <  1.0  and 

fh  +  3  A f  <  f  <  fh-  3A/ ,  where  f h  is  a  frequency  where  the  reference  rate  of  phase  change  is 
more  than  10°  /  Hz  and  A f  is  the  line  spacing  of  the  reference  spectra,  then  the  default  tolerance 
on  phase  will  be  +40°.  Otherwise  the  default  tolerance  on  phase  will  be  +10°. 

c.  Limiting:  See  paragraph  4.5.4. 

5.3  Laboratory  Data. 

In  the  case  the  reference  SDM  is  directly  employed  as  the  reference  in  a  MET  test  (i.e. 
input/output  (I/O))  Transformation  Control  as  discussed  in  reference  number  9),  it  is  not  a 
mathematical  based  hard  requirement  that  the  control  accelerometers  be  placed  in  the  exact  same 
location  in  the  laboratory  as  they  were  used  in  the  original  data  acquisition  phase.  The  critical 
parameter  is  that  all  control  locations  employed  in  the  laboratory  test  are  referenced  to  the  same 
“origin”  as  selected  in  the  original  VSD  development.  However,  it  is  often  desirable,  based  on 
making  position  specific  comparisons  between  field  and  laboratory  data,  to  match  the  laboratory 
control  locations  to  the  original  measurement  points. 
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6.  MDOF  VSD  METHODS. 

6.1  Options  Considered. 

Having  reviewed  the  data  acquisition  and  analysis  requirements,  this  section  is  dedicated  to 

17 

defining  the  steps  for  two  MDOF  VSD  methodologies  .  Method  I  is  processed  in  the  SDM 
domain  and  Method  II  conducts  averaging  steps  in  the  Cholesky  Domain.  An  example  follows 
in  paragraph  6.3. 

6.1.1  Method  I. 

The  following  is  a  10  step  outline  of  Method  I  (SDM  Domain)  MDOF  VSD: 

Step  1:  Prepare  to  convert  field  measurements  into  motion  DOFs. 

Identify  position  vectors  q  -r„  and  row  selection  vectors  e]  as  defined  in  paragraph 

4.4. 1.1,  corresponding  to  the  field  measurements. 

Identify  the  mission  scenario. 

Identify  the  frequency  bandwidth  of  interest. 

Identify  the  sampling  frequency  of  the  field  measurements. 

Step  2:  Transform  the  field  measurements  into  motion  DOF’s  per  equation  (4.3)  for  each  “Run” 
identified  in  the  mission  scenario. 

Step  3:  Compute  the  SDM  for  each  run  identified  in  Step  2.  The  dimension  of  the  resulting 
SDM’s  will  be  [6x6xd],  where  d  is  the  number  of  spectral  lines  being  considered, 
addresses  the  frequency  bandwidth  of  interest. 

Since  the  SDM  is  computed  from  measured  field  data,  it  should  be  positive  definite; 
however,  a  check  should  be  performed  to  verify  that  each  individual  SDM  is  positive 
definite.  This  serves  as  an  excellent  data  quality  check. 

Refer  to  the  guidance  in  Step  7  if  minor  corrections  are  required  to  force  an  individual 
SDM  to  be  positive  definite. 

Step  4:  Convert  the  CSD  terms  (the  off-diagonal  terms  of  the  SDM)  into  a  nonnalized  form  in 
which  the  magnitude  squared  of  the  cross  terms  correlates  to  the  ordinary  coherence 
while  leaving  the  phase  unchanged. 

This  is  accomplished  by  normalizing  (dividing)  the  CSD  terms  by  ^GXXGW  . 

While  it  is  not  absolutely  necessary  to  conduct  this  step,  it  is  often  easier  to 
understand  the  physical  meaning  of  the  CSD  terms  when  viewed  in  terms  of  phase 
and  coherence. 

Step  5:  Either  organize  all  of  the  SDM’s  for  the  Runs  of  interest  into  a  logical  structure  or  merge 
them  into  one  file  of  known  matrix  structure  such  as 

[SDM_Runl,SDM_Run2....SDM_RunN]  to  optimize  the  conduct  basic  statistics. 

Step  6:  Compute  a  weighted  average  SDM  of  the  N  SDM’s  of  Step  5. 
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It  is  critical  that  the  weighted  average  be  performed  on  the  true  complex  CSD  terms 
(not  the  nonnalized  SDM). 

The  weighting  factor  on  the  average  will  be  directly  correlated  to  the  mission  scenario 
times  identified  in  Step  1 .  If  the  individual  Runs  are  positive  definite,  the  resulting 
average  should  also  be  positive  definite.  However,  numerical  issues  may  yield  non¬ 
positive  definite  results.  To  minimize  numerical  issues,  average  only  the  lower 
triangular  portion  of  the  SDM  and  fill  in  the  upper  triangular  portion  of  the  SDM  by 
taking  advantage  of  the  Hermitian  structure  of  the  matrix  [16]. 

Any  type  of  enveloping  operation  should  be  avoided  as  it  is  highly  likely  to  yield  a 
non-positive  definite  result. 

Step  7:  As  SDM  data  are  manipulated  through  activities  such  as  averaging,  it  is  advisable  to 
verify  the  results  remain  positive  definite.  As  discussed  above,  occasional  numerical 
issues  may  be  of  concern  in  some  instances.  If  required,  force  the  SDM  computed  in 
Step  6  to  be  positive  definite. 

This  is  done  by  systematically  reducing  the  magnitude  of  the  cross  spectral  density 
tenns  until  the  Cholesky  decomposition  is  possible  at  each  depth  (spectral  line)  of  the 
SDM.  (If  required,  this  process  may  be  somewhat  conservative  in  its  reduction  of  the 
coherence  between  DOFs  in  that  the  systematic  reduction  of  cross  term  magnitudes  is 
applied  to  each  cross  term  equally). 

Step  8:  Scale  the  diagonal  tenns  of  the  autospectra  (the  diagonal  terms  of  the  SDM)  resulting 

from  Step  7  to  the  maximum  rms  level  of  each  of  the  N  SDM’s  in  Step  5  on  an  individual 
DOF  basis  using  Miner-Palmgren. 

Observe  that  a  new  total  test  time  will  be  computed  for  each  DOF  and  that  it  highly 
probably  that  the  resulting  test  times  for  each  DOF  will  not  be  the  same. 

Since  the  magnitude  of  the  autospectra  are  being  increased  while  not  modifying  the 
cross-spectral  density  terms,  the  resulting  scaled  SDM  should  still  be  positive 
definite.  However,  as  discussed  in  Step  7,  it  is  highly  recommended  that  anytime  a 
SDM  is  manipulated,  it  should  be  verified  that  the  resulting  SDM  remains  positive 
definite. 

Step  9:  Review  the  test  time  associated  with  each  DOF  resulting  from  Step  8  and  select  a 
reasonable  test  time  to  which  the  entire  SDM  may  be  referenced  to. 

Just  as  in  the  case  of  a  1-DOF  VSD  development,  one  should  consider  the  general 
guidance  to  keep  the  final  test  amplitudes  resulting  from  time  compression  to  be  no 
more  than  3  dB  above  the  maximum  measured  field  data.  Once  a  test  time  is 
selected,  reapply  Miner-Palmgren  as  required  per  DOF.  Again  make  sure  the 
resulting  SDM  is  positive  definite  and  modify  as  required  per  Step  7. 

Step  10:  Scale  the  results  from  Step  9  up  by  3  dB  to  account  for  uncontrolled  variables  such  as 
fleet  variations  and  scenario  conditions  not  considered  in  the  mission  scenario.  The 
resulting  SDM  and  the  test  time  association  per  Step  9  define  the  final  specification. 

This  is  accomplished  by  pre  and  post  multiplying  the  SDM  by  the  square  root  of  the 
ratio  of  the  desired  scaling  factor  as: 
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S Ynew  =  S,SYoMS,  (e-g- to  scale  the  SDM  ASD  terms  by  3  dB  while  keeping  the 
phase  and  ordinary  coherence  the  same,  the  diagonal  terms  of  Ss  would  be  defined  as 

ss,u=^y 

6.1.2  Method  II. 

The  following  is  a  10  step  outline  of  Method  II  (Cholesky  Domain)  MDOF  VSD: 

Steps  1-4:  Correlate  directly  to  Method  I  Outline. 

Step  5:  Perfonn  a  Cholesky  decomposition  on  the  individual  SDM  associated  with  each  Run  in 
the  mission  scenario. 

Since  each  individual  Run  was  based  on  a  physical  event,  the  individual  SDM’s 
should  be  positive  definite,  thereby  making  the  Cholesky  decomposition  possible. 
(Recall  all  Runs  would  have  been  tested  to  verify  each  was  positive  definite  or 
corrected  as  required  per  Step  3). 

Either  organize  all  of  the  lower  triangular  matrices  resulting  from  the  Cholesky 
decomposition  for  the  Runs  of  interest  into  a  logical  structure  or  merge  them  into  one 
file  of  known  matrix  structure  such  as  [CHOL_Runl,CHOL_Run2....CHOL_RunN] 
to  optimize  the  conduct  basic  statistics. 

Step  6:  Compute  a  weighted  average  Lower  Triangular  Matrix  of  the  N  Cholesky 
decompositions  of  Step  5. 

The  weighting  factor  on  the  average  will  be  directly  correlated  to  the  mission  scenario 
identified  in  Step  1 .  Note  that  the  resulting  average  will  still  consist  of  positive 
eigenvalues  implying  that  when  converted  back  into  the  SDM  fonnat  that  the  result 
will  be  positive  definite. 

Once  converted  back  into  the  SDM  domain,  the  resulting  CSD  tenns  will  generally  be 
highly  comparable  to  the  average  CSD  values  computed  in  Step  6  of  Method  I. 
However,  the  rms  levels  of  the  ASD  terms  will  not  be  the  same.  In  addition,  the 
spectral  shape  of  the  ASD  terms  will  generally  have  been  slightly  modified. 

Step  7:  Rescale  the  ASD  terms  of  the  SDM  resulting  from  Step  6  to  match  the  rms  levels  of  those 
in  Method  1  Step  6. 

Convert  the  CSD  tenns  (the  off-diagonal  tenns  of  the  SDM)  into  a  nonnalized  form 
in  which  the  magnitude  squared  of  the  cross  terms  conelates  to  the  ordinary 
coherence  while  leaving  the  phase  unchanged.  (Again,  while  it  is  not  absolutely 
necessary  to  conduct  this  step,  it  is  often  easier  to  understand  the  physical  meaning  of 
the  CSD  terms  when  viewed  in  tenns  of  phase  and  coherence). 

The  resulting  SDM  phase  and  coherence  are  expected  to  be  very  similar  to  those  of 
the  averaged  field  data  produced  in  Method  I.  The  ASD  terms  spectral  shapes  are 
expected  to  be  slightly  different  (i.e.  <  3  dB  per  spectral  line). 

Steps  8-10:  Correlate  directly  to  Method  I  Outline. 
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6.2  Example. 

a.  To  illustrate  the  process  discussed  above,  a  simple  example  was  derived  (Method  I  is 
addressed  first).  Using  an  available  wheeled  vehicle,  the  input  to  an  onboard  missile  storage 
rack  was  instrumented  as  shown  in  Figure  4.  The  transducer  at  the  center  of  Figure  4  was  placed 
at  the  user  defined  origin,  position  [0,0,0],  in  terms  of  a  Cartesian  coordinate  system.  In  a 
traditional  right  hand  orientation,  the  forward  direction  of  the  vehicle  was  defined  as  the  positive 
x-axis,  towards  the  vehicle  driver’s  side  was  considered  positive  y-axis,  and  through  the  vehicle 
roof  was  considered  the  positive  z-axis.  All  transducers  are  referenced  in  terms  of  their  relative 
placement  to  the  origin  as  discussed  previously  in  the  acceleration  transformation  section  of  this 
TOP. 


Figure  4.  Transducer  placement  (input  to  missile  rack). 


b.  Method  I  Example. 

(1)  Having  established  a  clear  coordinate  system  definition,  the  key  parameters 
discussed  in  Step  1  are  identified.  In  distance  units  of  inches,  the  positions  of  the  four  corner 
accelerometer  locations  used  in  this  example  are  defined  as: 

rl  =  [—17,— 6,0]',  r2  =  [-17,6,0]', r3  =  [17,  —6, 0] ' , r 4  =  [17, 6,0]'  ,  which  in  skew  symmetric  form  are: 
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For  convenience,  the  instrumentation  team  placed  the  tri-axial  transducers  such  that  the  channel 
used  to  measure  the  y-axis  motion  was  actually  180  degrees  out  of  phase  with  respect  to  the 
referenced  coordinate  system.  This  issue  is  addressed  by  simply  defining  row  selection  vectors 
ase*  =  [1,0,0],  =  [0,-10],  b  =  [0,0,1] .  Matrix  Ta  and  matrix  Ta  may  now  be  computed  as  per  the 


discussion  in  paragraph  4.4. 1.1  as: 


T  = 


0  0  0  0  6 

-1  0  0  0  17 

0  1  -6  17  0 

0  0  0  0  -6 

-1  0  0  0  17 

0  1  6  17  0 

0  0  0  0  6 

-1  0  0  0  -17 

0  1  -6  -17  0 

0  0  0  0  -6 

-1  0  0  0  -17 

0  1  6  -17  0 


T,  = 
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0 
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0 
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0.2500 
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0  0 
0.2500  0 
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-0.2500 

0 

0 

0 
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0 

0.2500 
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-0.0147 
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0.2500 

0 

0 

0 

0 

-0.0046 


0 

-0.2500 

0 

0 

0 

-0.0131 


0 

0 

0.2500 

0.0417 

-0.0147 

0 


The  field  data  were  sampled  at  4096  Hz  and  the  bandwidth  of  interest  is  500  Hz.  For  the 
example  at  hand,  a  mission  scenario  was  established  using  a  Beta  distribution  as  discussed  in 
reference  number  3,  and  is  illustrated  in  Table  1.  Allowing  for  the  time  associated  with  speeds 
below  5  miles  per  hour  (mph),  the  total  mileage  represented  is  approximately  300. 
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TABLET  MISSION  SCENARIO 


ROAD 

CLASSIFICATION 

SPEED 

(mph) 

TIME 

(hrs) 

DISTANCE 

(miles) 

Embedded  Rock 

5 

.690 

3.45 

10 

1.545 

15.45 

15 

.737 

11.05 

Cross  Country 

10 

5.18 

51.80  1 

20 

6.332 

126.64 

30 

2.002 

60.06 

Radial  Washboard 

5 

.811 

4.055 

7 

1.841 

12.88 

10 

1.183 

11.83 

(2)  The  field  data  were  then  converted  into  motion  DOFs,  {c}Motion ,  using 
Equation  4.3  per  Step  2. 

(3)  The  time  histories,  {c}Motion  were  then  transformed  into  the  frequency  domain  in 

the  form  of  a  SDM  per  run  as  described  in  Step  3.  Each  SDM  was  tested  per  the  Cholesky 
decomposition  property  and  verified  to  be  positive  definite. 

(4)  Each  SDM  was  then  normalized  as  suggested  in  Step  4  to  allow  the  analyst  to 
review  the  degree  of  coherence  between  DOFs. 

(5)  Per  Step  5,  the  SDMs  were  configured  into  a  convenient  structure  to  allow 
statistical  analysis.  MATLAB  was  employed  in  this  example  and  the  data  were  configured  as 
SDM_all=[SDM_Runl,SDM_Run2....SDM_Run8]  .  Observe  only  8  of  the  9  runs  identified  in 
the  scenario  are  being  considered.  In  reviewing  the  field  data,  the  5  mph  radial  washboard  data 
were  significantly  lower  than  the  rest  of  the  Runs,  detennined  to  have  no  effect  on  fatigue,  and 
were  not  considered  in  computing  the  basic  statistics  of  the  ensemble. 

(6)  Next,  per  Step  6,  a  weighted  average  in  terms  of  the  time  per  road  condition  as 
defined  in  Table  1  was  computed.  This  average  should  be  computed  in  terms  of  complex  CSD 
tenns,  not  the  nonnalized  SDM.  The  resulting  weighted  average  SDM  was  then  tested  at  each 
spectral  line  to  establish  whether  or  not  the  positive  definite  criterion  was  met.  Figure  5 
illustrates  the  weighted  average  SDM.  Taking  advantage  of  the  Hermitian  property  of  a  SDM, 
Figure  5  is  laid  out  such  that  the  lower  triangular  section  represents  the  phase  between  DOFs,  the 
upper  triangular  portion  represents  the  square  root  of  the  ordinary  coherence,  and  the  diagonal 
terms  are  the  ASDs  of  the  6  rigid  body  DOFs.  Although  too  small  to  review  in  detail  on  a  single 
page  as  shown,  the  coherence  plots  are  all  scaled  between  0.1  and  1.0.  This  is  to  illustrate  there 
is  some  level  of  coherence,  particularly  below  100  Hz  in  the  example  at  hand,  between  DOFs. 
Using  the  VSD  process  proposed,  the  analyst  will  try  to  keep  as  much  coherence  in  the  final 
specification  as  possible  while  still  ensuring  the  final  result  is  positive  definite. 
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Figure  5.  Normalized  weighted  average  SDM. 


(7)  In  order  to  address  the  possibility  of  having  to  deal  with  non-positive  definite 
results,  a  utility  was  written  which  gradually  and  equally  reduces  the  magnitudes  of  the  cross 
spectral  density  terms  until  the  positive  definite  criterion  is  met  per  Step  7.  This  technique 
actually  reduces  the  cross  term  magnitudes  of  some  CSDs  more  that  what  is  required. 

Addressing  this  potential  shortcoming  is  one  of  the  motivations  for  the  development  of 
Method  II. 

(8)  At  this  point,  per  Step  8,  the  rms  level  was  computed  for  each  ASD  (diagonal 
SDM  Entry)  over  the  bandwidth  of  interest  (3-500  Hz  in  this  example).  Each  ASD  was  then 
scaled  to  the  level  of  the  maximum  rms  level  via  Equation  4.6. 

(9)  Per  Step  9,  the  new  test  times  associated  with  each  ASD  were  also  documented. 

As  expected,  the  new  times  associated  with  each  DOF  were  no  longer  the  same.  Since  the  VSD 
effort  is  designed  to  yield  a  simultaneous  6-DOF  reference,  it  will  be  necessary  to  choose  a 
common  test  time  and  rescale  all  ASD  entries  to  the  selected  test  duration.  For  the  example  at 
hand,  a  test  duration  of  1 5  minutes  was  selected.  As  is  always  the  case  with  selection  of 
compressed  test  durations,  one  should  adhere  to  the  guidance  of  not  exaggerating  the  ASD  power 
levels  by  more  than  2:1.  Of  course  when  dealing  with  6  ASD  terms,  this  is  not  always  possible. 
In  such  cases,  the  analyst  should  avoid  increasing  the  dominant  DOFs  or  DOFs  with  known 
structural  shortcomings  by  more  than  3  dB  above  maximum  measured  ASD  levels. 
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(10)  The  terms  comprising  the  SDM  were  based  on  average  ASD  and  CSD  estimates, 
which  is  in  contrast  to  the  guidance  provided  in  reference  number  3,  in  which  the  ASD  levels 
carried  through  the  calculations  of  a  1-DOF  VSD  were  actually  based  on  an  ASD  computed  as 
the  sum  of  a  Mean  ASD  and  standard  deviation  computed  on  a  per  spectral  line  basis.  Working 
directly  with  the  mean  ASD  levels  is  intended  to  avoid  excessive  conservatism  in  the  VSD 
process.  Conservatism  intended  to  address  uncontrolled  variables  such  as  fleet  variations  and 
conditions  not  considered  in  the  mission  scenario  are  addressed  by  a  single  scalar  (+3  dB  in  this 
example)  in  Step  10.  Clearly  the  analyst  has  the  ability  to  modify  the  final  conservatism  level 
based  on  knowledge  of  the  specific  VSD  effort. 

The  final  reference  SDM  produced  by  Method  I  is  shown  in  Figure  6.  Observe  that  the  phase 
and  coherence  terms  are  essentially  unchanged  from  that  of  the  average  SDM  of  Figure  5. 


Figure  6.  Method  I  normalized  reference  SDM. 


c.  Method  II  Example.  The  first  four  steps  of  Method  II  correlate  directly  to  that  of 
Method  I.  The  major  deviation  in  Method  II  is  that  all  averaging  will  be  computed  in  the 
Cholesky  domain.  In  Step  5,  Cholesky  decompositions  are  carried  out  on  the  individual  SDM’s 
associated  with  each  Run  in  the  mission  scenario.  Since  each  individual  Run  was  based  on  a 
measured  physical  event,  the  individual  SDMs  were  positive  definite  as  expected,  thereby 
making  the  Cholesky  decomposition  possible.  In  the  event  that  a  given  Run  had  failed  the 
Cholesky  decomposition  and  all  measurement  locations  and  relative  polarities  were  verified; 
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investigate  the  spectral  lines  at  which  the  decomposition  fails.  If  the  decomposition  is  failing  at 
only  a  few  spectral  lines,  it  may  be  possible  to  salvage  the  measurement  employing  the  CSD 
magnitude  reduction  techniques  proposed  in  Method  I.  The  Cholesky  domain  data  were  then 
organized  into  a  convenient  structure  for  statistical  analysis.  As  in  Method  I,  Matlab  was  used  to 
compute  the  weighted  averages  and  the  Cholesky  domain  data  were  organized  as: 

CHOL_all=[CHOL_Runl,  CHOL_Run2 . ,CHOL_RunN]  .  In  Step  6,  a  weighted  average  in 

terms  of  the  time  per  road  condition  as  defined  in  Table  1  was  computed  over  the  lower 
triangular  matrix  of  the  eight  Cholesky  decompositions  of  Step  5.  The  weighted  average  was 
then  converted  back  into  the  SDM  domain.  As  expected,  the  coherence  characteristics  of  the 
resulting  SDM  were  comparable  with  that  of  Figure  5  and  the  rms  levels  of  the  ASD  terms 
required  rescaling  per  Step  7.  Steps  8-10  were  carried  out  directly  as  stated  in  the  Method  I 
outline. 

d.  The  reference  SDM  resulting  from  Method  II  (Figure  7)  yielded  similar  phase  and 
coherence  characteristics  to  that  of  the  reference  SDM  resulting  from  Method  I  (Figure  6).  Note 
that  the  Method  I  example  took  advantage  of  averaging  only  the  lower  triangular  CSD  terms, 
avoiding  potential  numerical  issues,  thereby  not  requiring  the  SDM  to  be  forced  positive  definite 
in  a  manner  that  would  result  in  lowering  the  coherence  in  a  more  conservative  manner  than 
required. 


Figure  7.  Method  II  reference  SDM. 
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e.  ASD  Comparisons.  Next,  the  minor  spectral  shape  deviations  between  the  ASD 
resulting  from  the  two  VSD  methods  discussed  will  be  illustrated.  Figures  8  and  9  show  the 
ASD  references  for  the  Z  axis  (vertical)  and  rotation  about  Z  axis  (Rz)  respectively,  as  produced 
from  both  VSD  methods.  The  ASD  references  are  superimposed  with  the  raw  (unexaggerated) 
reference  data  from  which  the  specifications  were  created.  Observe  that  the  ASD  shapes 
envelope  the  field  data  without  excessive  conservatism. 


Figure  8.  ASD  references  for  the  Z  axis. 
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Figure  9.  ASD  references  for  rotation  about  Z  axis  (Rz). 


f.  As  stated  previously,  the  test  duration  for  the  reference  SDM  yielded  by  both  Methods 
in  this  example  was  established  to  be  15  minutes.  Clearly,  as  illustrated  in  Figures  8  and  9  the 
associated  ASD  references  are  highly  correlated. 

6.3  Concluding  Remarks. 

a.  Two  techniques  were  defined  for  establishing  an  input  specification  for  a  MDOF 
system.  It  was  shown  that  simple  enveloping  techniques  are  not  appropriate  when  considering 
CSD  terms  due  to  the  sensitivity  of  such  operations  associated  with  maintaining  a  physically 
realizable  reference.  The  resulting  SDM  references  yielded  through  the  process  outlined  are 
fully  populated  SDM’s.  Importing  the  fully  populated  SDM  into  the  MDOF  control  system  in  an 
efficient  manner  is  essential  due  to  the  volume  of  information  involved. 

b.  While  synthesizing  a  drive  signal  with  CSD  characteristics  of  the  field  data  is  desired, 
it  is  recognized  that  the  mechanical  impedance  of  the  laboratory  configuration  is  highly  unlikely 
to  match  that  of  the  field  data.  Therefore,  it  will  be  difficult  to  maintain  CSD  characteristics 
across  the  spectral  bandwidth  of  interest  and  thus,  the  control  hierarchy  will  generally  place 
emphasis  on  the  ASD  terms.  Also,  it  is  not  uncommon  in  MDOF  tests  for  a  specific  mechanical 
degree-of-freedom  to  consist  of  a  very  small  percentage  of  the  composite  energy  across  all 
mechanical  degrees-of-freedom.  In  such  cases,  the  associated  error  for  the  low  DOF  will  often 
be  higher  than  the  desired  test  tolerances  and  considering  global  test  tolerances  may  need  to  be 
considered. 
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c.  Care  was  taken  in  the  examples  provided  to  limit  the  amount  of  conservatism  in  the 
VSD  process.  One  quickly  realizes  that  the  amount  of  conservatism  is  cumulative  across  degrees 
of  freedom  and  if  not  managed  carefully  will  yield  test  levels  significantly  higher  than  the 
measured  environment.  Unlike,  the  common  technique  of  essentially  adding  3  dB  to  all 
measurements  prior  to  conducting  averaging  or  enveloping  techniques  in  the  1-DOF  arena  per 
reference  number  3,  all  weighted  averages  in  the  6-DOF  examples  shown  were  based  on  raw 
averaged  data.  Conservatism  to  account  for  variables  such  as  fleet  variability  and  mission 
scenario  omissions  were  added  in  the  final  step.  Magnitude  amplification  associated  with  time 
compression  techniques  was  limited  to  no  more  than  maximum  measured  levels.  Also,  on  the 
subject  of  tolerances,  one  may  find  it  reasonable  to  define  phase  and  coherence  tolerances  over 
only  a  portion  of  the  test  bandwidth.  In  the  example  provided,  the  coherence  dropped  off 
considerably  at  frequencies  above  50  Hz.  Since  the  phase  term  is  essentially  a  random  variable 
for  low  coherence,  setting  tolerances  for  frequencies  greater  than  50  Hz  would  not  be 
recommended  for  the  example  shown. 
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Term 


Definition 


Laboratory 
Vibration  Test 
Schedule  (LVTS) 


Multi- Axis  (MA) 


Multi-Degree-of- 
Freedom  (MDOF) 


Multiple 

Exciter/Multiple- 
Axis  (MEMA) 


Multiple-Exciter/ 
Single- Axis 
(MESA) 


Multiple- 
Input/Multiple- 
Output  (MIMO) 


All  information  required  to  perform  a  vibration  test  on  a  vibration 
exciter.  Information  typically  includes:  a  broadband  spectra  (or 
profile),  sine  or  narrowband  infonnation  (if  used),  test  run  time, 
control  accelerometer  locations,  control  methods  and  tolerances, 
and  any  test  specific  infonnation  required. 

Excitation  or  response  measurement  that  requires  more  than  one 
unique  vector  for  description. 

Motion  defined  by  test  item  movement  along  or  about  more  than 
one  axis  whose  description  requires  two  or  more  coordinates  to 
completely  define  the  position  of  the  item  at  any  instant. 

Application  of  multiple  exciters  providing  dynamic  input  to  the  test 
item  in  a  way  that  requires  more  than  a  single  vector  for  complete 
description  of  excitation  and  measurement.  Note  that  many  multi¬ 
axis  test  platfonn  configurations  have  been  built  in  recent  years. 
Many  are  optimally  configured  such  that  the  number  of  mechanical 
DOFs  is  equivalent  to  the  number  of  exciters.  There  are  also  over¬ 
actuated  systems  consisting  of  more  exciters  than  mechanical 
DOFs.  In  each  case,  the  dynamic  properties  vary  between  designs 
and  must  be  considered  in  the  design  of  a  MET. 

Application  of  multiple  exciters  providing  dynamic  input  to  the  test 
item  in  a  single  vector  direction.  For  example,  extended  materiel 
might  require  excitation  at  the  forward  and  aft  end  in  a  single  vector 
axis.  For  the  case  in  which  the  two  exciters  are  driven  to  a  common 
specification  with  respect  to  both  phase  and  amplitude,  the  output 
may  be  described  basically  in  the  one  axis  of  excitation.  For  the 
case  in  which  the  two  exciters  are  driven  to  independent  magnitude 
and/or  phase  specifications,  the  output  may  need  to  be  described  in 
terms  of  a  forward  axis  and  aft  axis  and,  perhaps,  a  rotational  axis 
about  the  test  item’s  center-of-gravity  (CG). 

Refers  to  input  of  a  multiple  drive  signals  to  an  exciter  system 
configuration  in  a  MDOF  configuration,  and  multiple  measured 
outputs  from  the  fixture  or  test  item  in  a  MDOF  configuration.  It  is 
important  to  note  that  generally  there  is  no  one-to-one 
correspondence  between  inputs  and  outputs,  and  the  number  of 
inputs  and  number  of  outputs  may  be  different. 
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Term 

Multiple- 
Input/S  ingle- 
Output  (MISO) 


Scenario 


Single- Axis  (SA) 


Single-Degree-of- 
Freedom  (SDOF) 


Single- 

Exciter/Single- 
Axis  (SESA) 

Single- 

Input/Multiple- 
Output  (SIMO) 


Single- 
Input/S  ingle- 
Output  (SISO) 


Definition 

Refers  to  input  of  a  multiple  drive  signals  to  an  exciter  system 
configuration  in  a  MDOF  configuration  and  a  single  measured 
output  from  the  fixture  or  test  item  in  a  SDOF  configuration.  This 
tenninology  is  most  used  in  measurement  data  processing  where  the 
single  output  is  a  composite  of  measurements  from  multiple  inputs. 

A  tabulation  of  expected  exposure  events  and  the  corresponding 
durations. 

Excitation  or  response  measurement  in  a  unique  single  vector 
direction  (linear  or  rotational).  For  rotational  axis,  the  vector 
direction  is  perpendicular  to  the  plane  of  rotation  of  the  exciter  or 
test  item. 

Motion  defined  by  materiel  movement  along  or  about  a  single  axis 
whose  description  requires  only  one  coordinate  to  completely 
define  the  position  of  the  item  at  any  instant. 

Application  of  a  single  exciter  providing  dynamic  input  to  the  test 
item  in  a  single  vector  direction. 


Refers  to  input  of  a  single  drive  signal  to  an  exciter  system  in  a 
SDOF  configuration,  and  multiple  measured  outputs  from  the 
fixture  or  test  item  in  a  MDOF  configuration.  In  general,  for 
specification  purposes  the  dynamic  behavior  of  the  test  item  will 
not  be  assumed  to  contribute  to  the  output  DOF,  i.e.,  measured 
rotation  of  an  extended  test  item  that  is  being  excited  in  a  cantilever 
mode  will  still  basically  be  considered  as  a  SET  with  linear 
acceleration  characterizing  the  output. 

Refers  to  input  of  a  single  drive  signal  to  an  exciter  system  in  a 
SDOF  configuration  and  a  single  measured  output  from  the  fixture 
or  test  item  in  a  SDOF  configuration. 
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AECTP 

ASD 


CG 

CSD 

dB 

DFT 

DOF 

DTC 

FRF 

g/v 

Hz 

I/O 

IEEE 

IES 

IEST 

LCEP 

LOA 

LVTS 

MA 

MDOF 

MEMA 

MESA 

MET 

MIL-STD 

MIMO 

MISO 

NATO 

NBROR 

pinv 

PSD 

rms 

RTC 
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Allied  Environmental  Conditions  Test  Publication 

auto  spectral  density  (also  referred  to  as  the  power  spectral  density 

(PSD) 

center  of  gravity 
cross  spectral  density 

decibel 

discrete  Fourier  transform 
degree  of  freedom 

US  Army  Developmental  Test  Command 
frequency  response  function 
gravitational  units/volts  of  drive 
hertz 

input/output 

Institute  of  Electrical  and  Electronics  Engineers 

Institute  of  Environmental  Sciences 

Institute  of  Environmental  Sciences  and  Technology 

Life  Cycle  Environment  Profile 
line  of  action 

Laboratory  Vibration  Test  Schedule 
multi-axis 

multiple  degree-of-freedom 
multiple-exciter  multiple-axis 
multiple-exciter  single-axis 
multiple  exciter  test 
Military  Standard 
multiple-input  multiple-output 
multiple-input  single-output 

North  Atlantic  Treaty  Organization 
narrowband  random  on  random 

Moore  Penrose  pseudo  inverse 
power  spectral  density 

root  mean  square 

US  Army  Redstone  Test  Center 
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SA 

single-axis 

SDM 

spectral  density  matrix 

SDOF 

single  degree-of-freedom 

SESA 

single-exciter/ single-axis 

SIMO 

single-input  multiple-output 

SISO 

single-input  single-output 

SOR 

sine-on-random 

TOP 

Test  Operations  Procedure 

TWR 

Time  Waveform  Replication 

UUT 

unit  under  test 

VSD 

Vibration  Schedule  Development 
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APPENDIX  C.  NOMENCLATURE. 

Definition 

A  vector  where  each  element  is  a  discrete  time  history  or  function 
of  frequency,  the  discrete  Fourier  transfonn  (DFT)  of  a  time 
history.  In  general  lower  case  letters  will  be  used  for  functions  of 
time  and  upper  case  letters  will  be  used  for  functions  of  frequency. 
Sometimes  lower  case  letters  are  used  to  designate  an  element  in  a 
matrix. 

Will  denote  a  matrix.  Usually  a  third  dimension  will  denote  time 
samples  or  samples  as  a  function  of  frequency. 

The  complex  conjugate  transpose  of  a  matrix. 

The  Moore  Penrose  pseudo  inverse  of  a  matrix. 

The  transpose  of  a  matrix. 

Over  a  variable  will  denote  an  estimated  value. 

The  vector  of  return  acceleration  signals. 

The  spectral  density  matrix  of  the  return  signals,  typically  in  units 
of  g2 / Hz. 

A  vector  of  the  control  signals  from  a  MIMO  system.  Each  element 
in  the  vector  is  a  function  of  time.  It  can  be  thought  of  as  a  2 
dimensional  matrix.  First  dimension  is  the  input  number.  The 
second  dimension  is  the  time  index. 

The  DFT  of  {c}. 

The  spectral  density  matrix  of  the  control  signals.  The  diagonal 
elements  are  the  real  auto-spectral  densities  of  the  control  signals. 
The  off  diagonal  elements  are  complex  functions  of  frequency 
giving  the  cross  spectral  density  between  pairs  of  control  signals. 

A  vector  of  drive  signals  into  a  MIMO  system.  Each  element  in  the 
vector  is  a  function  of  time.  It  can  be  thought  of  as  a  2  dimensional 
matrix.  First  dimension  is  the  input  number.  The  second 
dimension  is  the  time  index. 
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Term 

[D] 

E[] 

g 

[H] 


Ns 

Nc 

Na 


Nd 


R 

W 

S 

ScD 


Definition 

The  drive  signals  in  the  frequency  domain.  \d\  is  formed  from  [D] 
using  a  method  called  time  domain  randomization.  Initially 
D  =  ZRZ'. 

The  expected  value. 

The  acceleration  of  gravity. 

A  matrix  of  frequency  response  functions  (FRF’s)  relating  the 
control  system  response  to  the  drive  signals.  Typically  the  elements 
will  have  units  of  g/V.  Each  element  is  a  frequency  response 
function.  A  third  dimension  typically  is  the  amplitude  as  a  function 
of  a  set  of  frequencies  relating  to  the  DFT  of  the  input  and  response 
signals. 

The  number  of  drive  signals,  the  number  of  shakers. 

The  number  of  control  signals. 

The  number  of  acceleration  return  signals. 

The  number  of  output  control  variables. 

The  reference  control  spectral  density  matrix;  the  desired  spectral 
density  matrix. 

The  vector  of  shaker  drive  voltages. 

The  spectral  density  matrix  of  the  drives  in  shaker  space. 

The  spectral  density  matrix  between  the  control  signal  and  the 
drives  to  the  shakers. 


Ta  The  acceleration  to  control  space  transformation  matrix. 

Ts  The  drive  in  the  control  space  to  voltages  {s}  to  the  shakers 

transformation  matrix. 

Z=  H'  The  system  impedance  matrix,  typically  units  of  volts/g. 
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APPENDIX  D.  MATRIX  ALGEBRA  REVIEW. 

A  matrix  is  an  array  of  numbers  arraigned  in  rows  and  columns.  The  size  of  the  matrix  is 
typically  stated  an  [n,m]  or  n  x  m,  where  n  is  the  number  of  rows  and  m  is  the  number  of 
columns.  In  this  document  3  dimensional  matrices  are  also  used  where  the  third 
dimension  is  typically  samples  in  either  the  time  or  frequency  domain.  This  Appendix 
will  discuss  only  two  dimensional  matrices.  It  is  assumed  that  if  the  matrix  has  3 
dimensions,  that  the  operations  can  be  perfonned  on  each  2  dimensional  matrix  along  the 
third  dimension.  For  example  if  the  matrix  is  a  matrix  of  frequency  response  functions, 
matrix  operations  will  be  perfonned  at  each  frequency  line.  The  definitions  provided  in 
this  appendix  are  based  on  information  provided  primarily  in  reference  numbers  12  and 
13. 


a.  SDM:  A  spectral  density  matrix  is  a  3  dimensional  matrix.  At  each  frequency 
line  (the  3rd  index)  the  matrix  is  a  square  complex  matrix.  Each  diagonal  element  is  the 
autospectrum  of  the  corresponding  element.  Loosely  an  element  in  the  SDM  is  defined 
as: 


Gn  (k)  =  2  lim ±E[X .  (k,  T)X]  (k,  T)] 

J  T— »oo  y  J 


where:  G  Jk)  is  the  cross  spectral  density  between  the  j’th  and  i’th  random  processes. 
X  j(k,T)  and  A  (k,  T)  are  the  discrete  Fourier  transforms  of  the  time  histories,  and  k  is 

the  frequency  index.  If  i  =  j,  the  spectrum  is  called  the  autospectrum  or  the  power 
spectrum.  In  reality,  the  true  spectral  density  is  generally  not  known  and  an  estimate  is 
employed.  Some  authors  define  the  elements  as: 

Gy  (k)  =  2  lim  i  E[  x;  (k,  T)X  j  (k,T)\ 

J  r— >oo  y  j 


The  SDM  matrix  is  Hermitian  positive  definite. 

b.  Hermitian  Matrix:  A  matrix,  A,  is  Hermitian  if  the  diagonal  elements  are  real 
positive  numbers  and  the  corresponding  off  diagonal  elements  are  complex  conjugate 
pairs: 


an  =  positive  real  number 
a  ^  =  a*.  =  conjiUj) 

where:  a  ..  is  the  element  fonn  j’th  row,  i’th  column  of  A. 
Note:  All  valid  spectral  density  matrices  (SDM)  are  Hennitian. 
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c.  Positive  Definite  Matrix  and  Positive  Semi-Definite  Matrix:  If  a  square 
Hennitian  matrix,  A,  has  all  positive  eigenvalues,  the  matrix  is  positive  definite.  If  the 
matrix  has  zero  eigenvalues  the  matrix  is  positive  semi-definite.  A  Cholesky 
decomposition  is  possible  for  all  positive  definite  matrices. 

A  =  LL' 

where:  L  is  a  lower  triangular  matrix  with  real  positive  values  on  the  diagonal.  L’  is  the 
complex  conjugate  transpose  of  L.  If  the  matrix,  A,  is  positive  semi-definite,  special  care 
must  be  taken  in  computing  L.  If  a  zero  element  is  found  on  the  diagonal  of  L,  the  entire 
column  must  be  set  to  zero.  Computing  the  Cholesky  decomposition  is  actually  the 
easiest  way  to  check  for  positive  definite.  If  the  algorithm  fails  the  matrix,  A  is  not 
positive  definite. 

d.  Transfonnation  of  a  Positive  Definite  Matrix: 

Let  B  =  HAH' 

If  the  matrix  A  is  positive  definite,  B  is  positive  definite. 

Note:  All  valid  SDMs  are  positive  semi-definite  or  positive  definite.  Because  some 
noise  is  always  present  in  measured  data,  a  measured  SDM  will  always  be  positive 
definite. 


2 

e.  Ordinary  Coherence,  y  :  The  ordinary  coherence  between  two  signals  is 
defined  as: 


rh  = 


'12 


GnG 


2 


22 


Gn  is  the  cross  spectral  density  between  the  signals  and  Gn  and  G22  are  the  two 
autospectra. 

The  ordinary  coherence  is  bounded  by  0  <  rh  ^  i  • 

Coherence  is  a  measure  of  the  linear  relationship  between  the  signals.  If  the  coherence  is 
unity,  a  perfect  linear  relationship  exists  between  the  signals.  If  the  coherence  is  zero,  the 
signals  are  said  to  be  independent,  and  there  is  no  linear  relationship  between  the  signals. 

If  one  or  more  of  the  ordinary  coherences  in  a  SDM  are  in  unity  at  any  frequency,  the 
matrix  is  positive  semi-definite  at  that  frequency. 
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f.  Singular  Value  Decomposition:  Singular  value  decomposition  has  several 
applications  in  MIMO  testing.  Singular  value  decomposition  is  defined  as: 

M=USV’ 

M  is  any  matrix.  U  and  V’  are  orthonormal.  This  implies  that: 


UU’=I  and  VV’=I 


S  is  a  diagonal  matrix  of  non-negative  real  numbers.  A  common  convention  is  to  order 
the  diagonal  elements  of  S  in  a  non-increasing  fashion. 

g.  Pseudo  inverse:  The  Moore  Penrose  pseudo  inverse  is  used  often  in  MIMO 
control.  Some  of  the  properties  are  discussed  below.  The  Moore  Penrose  pseudo  inverse 
can  be  derived  as  follows: 


M  =  USV' 

U'M  =  U'USV'=  SV' 

S_1U'M  =  S_1SV’  =  V’ 

VS'U'M  =  W'=  I 

Mf  =  VS_1U'  is  known  as  the  pseudo  inverse  of  M. 

The  inverse  of  the  reduced  S  is  easy  since  the  matrix  is  diagonal.  To  compute  S'1  the 
elements  greater  than  a  tolerance  are  inverted  and  kept,  the  elements  less  than  a  tolerance 
are  replaced  by  zero. 

MMf  M  M  and  MfM  Mf  =  Mf 
MM1  and  M^M  are  Hermitian 

If  the  number  of  columns  in  M  exceed  the  number  of  rows  and  the  rows  are  independent 
MM-I.  If  the  number  of  rows  in  M  exceeds  the  number  of  columns  and  the  columns  are 
independent  M^M=I.  For  a  more  complete  discussion  see  the  help  file  for  pinv  in 
MATLAB. 

h.  Matrix  Rank:  The  rank  of  a  matrix,  M,  equals  the  number  of  non-zero  singular 
values  in  M.  In  numerical  linear  algebra,  the  singular  values  can  be  used  to  determine  the 
effective  rank  of  a  matrix.  Define  a  measure  of  singular  values  as  the  ratio  of  the  singular 
values  and  the  largest  singular  value.  Let  r  be  the  number  values  greater  than  a  threshold. 
Where  the  measure  is  less  than  the  threshold,  set  the  singular  values  to  zero.  The  number 
of  non-zero  singular  values  in  the  resulting  matrix  is  the  effective  rank  of  the  matrix.  The 
effective  rank  of  the  matrix  is  r.  For  a  square  matrix,  if  r  is  less  than  the  number  of  rows 
and  columns  in  the  matrix,  the  matrix  is  said  to  be  ill  conditioned. 
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i.  Matrix  Approximation:  Let  M  =  usv' 

where:  s  =  a  diagonal  matrix  of  the  singular  values  greater  than  a  threshold  defined  as  the 
ratio  of  the  singular  values  divided  by  the  largest  singular  value.  Let  n  =  the  number  of 
kept  singular  values,  s  has  n  rows  and  columns,  u  is  the  first  n  columns  of  U.  v’ is  the 
first  n  rows  of  V’. 


M  minimizes  LS-wsvl 


Hence,  M  is  a  very  useful  approximation  of  M. 

j.  Frobenius  Norm:  The  Frobenius  Norm  of  matrix  M  is  defined  as: 


M,=,E  X  mu  =y/trace(A'A)= 


min  (m,n) 


\  i=\  j= 1  V  i= 1 

where:  cr  are  the  singular  values  of  M. 

k.  Trace:  The  trace  of  a  positive  definite  matrix  is  defined  as  the  sum  of  the 
diagonal  elements.  An  important  property  of  the  trace  often  of  use  is: 

trace(AB)  =  trace(BA) 


1.  Rescaling  the  Autospectra:  When  generating  a  SDM  it  might  sometimes  be 
useful  to  rescale  the  autospectra  and  be  assured  that  the  result  remains  positive  definite. 
This  can  be  accomplished  by  pre  and  post  multiplying  by  a  diagonal  matrix  of  scaling 
factors.  The  triple  product  will  rescale  the  autospectra  while  keeping  the  coherence  and 
phase  between  pairs  of  channels  unchanged. 


Gnew=SGoldS' 


where:  Gnew  is  the  new  positive  definite  SDM,  Gold  is  the  original  positive  definite 
SDM,  and  S  is  a  diagonal  matrix  of  scaling  factors.  Each  autospectra  will  be  scaled  by 
the  corresponding  element  in  S2 . 

This  is  a  convenient  way  to  generate  the  nonnalized  SDM  (the  diagonal  elements  are  the 
autospectra  and  the  magnitude  squared  of  the  off  diagonal  tenns  are  the  ordinary 
coherence  and  the  phase  is  the  phase  of  the  cross  spectra).  The  normalized  form  is 
computed  by  rescaling  the  SDM  to  unity  autospectra  by  pre  and  post  multiplying  the 
SDM  by  a  diagonal  matrix  whose  terms  are  the  inverse  square  root  of  the  autospectra. 
The  resulting  unity  autospectra  are  then  replaced  by  the  original  autospectra. 
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The  inverse  is  computed  by  replacing  the  diagonal  autospectra  by  ones  and  then  rescaling 
by  pre  and  post  multiplying  by  a  diagonal  matrix  whose  terms  are  the  square  root  of  the 
original  autospectra. 

m.  Proof  1 : 

An  element  in  Sy  is  given  by,  where  n  =  number  of  inputs,  and  m  =  number  of  outputs 

i  =  \:m  j  =  l:m 

r= 1  k= 1 

A  diagonal  element  is  given  by: 


r„ = X’£x„hi,h;  =  YLX*  i  H>  f  = 1  ■ m 


r= 1  k= 1 


r= 1  k= 1 


If  Sx  is  diagonal,  =  0,  if  r  ^  k,  (a-3)  reduces  to: 


i  =  1 :  /n 


This  can  be  written  as  a  set  of  linear  equations: 

Y  =  HX 

Which  can  be  solved  for  X  as:  X  =  H  1 Y 
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